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In a method for isolating anatomical structures contained within a three-dimensional data set. a morphological skeleton (302) is first 
fomied from the three-dimensional data set. Next, a seed data point (310) is selected from within the morphological skeleton. The seed 
data point is contained within a desired anatomical structure to be displayed and/or analyzed. Fuzzy connectivity (314) is utilized to define 
additional data points of the desired anatomical structure, so as to facilitate reconstruction of substantially only the desired anatomical 
structure. Such reconstruction of substantially only the desired anatomical structure facilitates viewing and analysis thereof by reducing 
complexity of the image and eliminating obstructing tissue. 
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DATA DIMENSIONAL SIEVING AND FUZZY CONNECTIVITY 
FOR MRI IMAGE ANALYSIS 

Field of t:he Inveniiion 

The present invention relates generally to medical 
imaging systems and digital signal processing. It relates 
more particularly to the use of data dimensional sieving 
and fuzzy connectivity to facilitate analysis and review of 
three-dimensional medical images such as those produced by 
magnetic resonance imaging (MRI) devices and the like- 

Background of the Invention 
Tomographic imaging techniques for use in medical 
applications are well known. Examples of such techniques 
include magnetic resonance imaging (MRI), computer aided 
tomography (CAT), and positron emission tomography (PET). 
In each of these techniques, a plurality of cross-sectional 
two-dimensional images, i.e., slices, of a body portion are 
generated and processed so as to provide a three- 
dimensional model of the imaged body portion. 

To utilize such a three-dimensional model of a body 
portion, slices or images taken along planes of interest 
are generated and then images are printed or otherwise 
displayed for viewing. The slices viewed may be at various 
angles with respect to the three-dimensional model. They 
do not necessarily correspond to the angles of the slides 
from which the three-dimensional model was originally 
constructed. Thus, the medical diagnostition may review 
images taken along any desired plane within the three- 
dimensional model. This provides for a great deal of 
flexibility in the use of the three-dimensional model as a 
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diagnostic tool. 

As those skilled in the art will appreciate, such 
three-dimensional images provide a valuable tool to the 
medical diagnostition in a manner which is non-invasive, 
5 and which is therefore considered to be of very low risk to 
the patient. 

However, although such three-dimensional imaging 
techniques have proven extremely useful for their intended 
purposes, they still possess inherent deficiencies which 

10 detract from their overall effectiveness. More 
particularly, it is frequently difficult to interpret the 
viewed two-dimensional slices or images when the anatomical 
structures of interest are surrounded by and/or intermixed 
with various other anatomical structures. The undesirable 

15 presence of such superfluous imagery only complicates the 
image, making it much more difficult to view and interpret 

the desired imagery. 

For example, viewing delicate portions of the vascular 
system is typically difficult since veins, arteries, and 
20 capillaries are intermixed with surrounding tissue. This 
makes it very difficult to distinguish the desired portions 
of the vascular system from surrounding tissue. Often, 
only slight changes in the intensity of the image 
distinguish a desired anatomical structure from surrounding 

25 tissue. , * 

Thus, it is desirable to provide a method for 
isolating anatomical structures of interest such that 
surrounding tissue is not displayed along therewith. In 
this manner, the medical diagnostition may view only the 

30 unobstructed anatomical structures of interest. This 
vastly reduces the complexity of the image and thus 
minimizes confusion as to precisely what portions of the 
image relate to the anatomical structure of interest. 



35 
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Summary of t:he Invent:ion 
In the analysis and review of three-dimensional 
medical imaging, it is of critical importance to be able to 
measure and analyze image features having various fractal 
dimensionalities from zero dimensions to three dimensions. 
For example, veins and arteries are characterized as one- 
dimensional curvilinear forms, while capillaries exhibit 
one plus fractal dimensions, typically exhibiting 
fractional fractal dimensionality. Tumors have three- 
dimensional fractal forms and exhibit smaller fractal 
dimensions if metastases are considered. 

The present invention specifically addresses and 
alleviates the above-mentioned deficiencies with the prior 
art. More particularly, the present invention comprises a 
method for isolating anatomical structures contained within 
a three-dimensional data set, e.g., a three-dimensional 
model formed by MRI, a CAT scan, or a PET scan. The method 
comprises the steps of forming a morphological skeleton of 
the three-dimensional data set, selecting a seed data point 
within the morphological skeleton so as to identify a 
desired anatomical structure to be displayed or analyzed, 
and utilizing fuzzy connectivity to define additional data 
points of the desired anatomical structure so as to 
reconstruct substantially only the desired anatomical 
structure. Reconstruction of substantially only the 
desired anatomical structure facilitates the review and 
analysis of the anatomical structure. 

For example, if it is desirable to obtain a three- 
dimensional data set containing only data points which are 
representative of the brain, then the patient's head may be 
imaged via MRI, CAT, PET scanning techniques or the like to 
provide a three-dimensional model of substantially the 
entire head. The three-dimensional data set which defines 
this model is then processed so as to form a morphological 
skeleton thereof. 



wo 97/09690 



PCT/US96/14500 



-4- 

An operator then selects a seed data point within the 
morphological skeleton corresponding to the patient's 
brain. This is typically accomplished by viewing the 
morphological skeleton on a display such as a CRT. The 
morphological skeleton maintains all of the data available 
in the original three-dimensional data set. However, in 
the morphological skeleton, anatomical structiires are 
separated from one another, based upon the fractal 
dimensionality thereof . Thus , anatomical structures having 
a fractal dimensionality of less than one dimension are 
separated from those having a fractal dimensionality of 
less than two dimensions and the anatomical structures are 
separated from those having a fractal dimensionality of 
less than three dimensions. 

After selecting a seed data point within the brain, 
fuzzy connectivity is utilized to define the additional 
data points which are required to provide a substantially 
complete image of the brain. Reconstruction of the brain 
is simply the reverse of the process utilized to form the 
morphological skeleton. With the use of fuzzy connectivity 
to define the set of points defining the brain, it appears 
that all of the features thereof are substantially utilized 
in the reconstruction process. Reconstruction of the brain 
without the use of fuzzy connectivity would result in the 
loss of substantial surface details thereof. For example, 
the surface texture and even, to a lesser degree, the 
convolutions of the brain, would tend to be degraded or 
smoothed . 

The morphological skeleton is formed by recursive 
opening and erosion of the three-dimensional data set so as 
to form a plurality of residuals which define the 
morphological skeleton. Reconstructing a desired anatomic 
structure from the morphological skeleton comprises 
performing the opposite procedure from that utilized to 
form the morphological skeleton. Thus, reconstruction 
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comprises recursive dilation and closing of the 
morphological skeleton. As those skilled in the art are 
aware, each step of opening comprises an erosion followed 
by a dilation and each step of the closing comprises a 
5 dilation followed by an erosion. 

The use of fuzzy connectivity during the 
reconstruction process assures that substantially all of 
the data points associated with the desired anatomical 
structure are utilized in the reconstruction process. 
0 According to the preferred embodiment of the present 

invention, a generally spherical structuring element is 
utilized in both the formation of the morphological 
skeleton and the reconstruction process. However, those 
skilled in the art will appreciate that various other 
.5 shapes of structuring elements are likewise suitable. 
Indeed, it has been found that various different shapes of 
structuring elements are particularly suited for use with 
various different dimensionalities or shapes of anatomical 
structures . 

0 Generally, the seed data point is selected by 

positioning a cursor at a desired point on an image being 
displayed upon a monitor. Thus, the operator may simply 
visually identify and manually select a seed within the 
organ or anatomical structure of interest. However, as 

!5 those skilled in the art will appreciate, various different 
computer algorithms may be utilized in the selection of 
such a seed. For example, the operator may simply initiate 
an algorithm which selects the largest organ within a given 
volume. Thus, if the operator desires to select the brain 

30 for reconstruction, the operator could merely select the 
largest organ within the head. 

The use of fuzzy connectivity to define additional 
data points of the desired anatomical structure comprises 
defining connectivity based upon the size and shape of a 

35 structuring element utilizing a fuzzy generalization of 
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mathematically defined distances between sets of data 
points as a criterion. This is accomplished based upon a 
modified Hausdorff metric. 

Thus^ separation of such anatomical features from one 
another according to the present invention is accomplished 
via dimensional sieving. Dimensional sieving results in 
the formation of a morphological skeleton utilizing the 
recursive opening and erosion processing according to well 
known principles. The opening and erosion processes are 
described in detail in "Morphological Systems for Multi- 
Dimensional Signal Processing" by Petros Maragos and Ronald 
W- Schafer, Proceeds of the IEEE > Volume 78, No. 4, April 
1990; "Morphological Filters-Part I: Their Set-Theoretic 
Analysis and Relation to Linear Shift-Invariant Filters", 
by Petros Maragos and Ronald W. Schafer, IEEE Transactions 
on Acoustics. Speech, and Sign al Processing, Volume ASSP- 
36, No. 8, August 1987; and "Morphological Filters, Part 
II: Their Relations to Median Order-Statistic, and Stack 
Filters", by Petros Maragos and Ronald W. Schafer, IEEE 
Transactions on Acoustics. Sp eech, and Signal Processing, 
Volume ASSP-35, No. 8, August 1987. 

According to the present invention, a cascade of data 
dimensional sieving filters are used directly with a three- 
dimensional image from an MRI device or the like to isolate 
structures such as arteries and veins from surrounding 
tissue for unobstructed visualization. This cascade of 
data dimensional sieving filters comprises the use of a 
generally spherical structuring element, followed by the 
use of a two-dimensional surface structuring element, 
followed by the use of a curvilinear structuring element, 
followed by the use of a point structuring element. 

Thus, to provide for the identification of desired 
dimensional features within the multi-dimensional data set 
provided by a tomographic imaging device, a data 
dimensional sieving algorithm separates the data based upon 
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the dimensional characteristics 
structures contained therein. ^ 



5 of the anatomical 
The algorithm utilizes 



filters which resemble geometric constructions such as 
lines, disks, and spheres, to sieve multi-dimensional 
features of curves, surfaces, and regions, as well as 
features of fractal dimensions in between. 

A hierarchy of dimensional filters is thus utilized to 
first remove features of less than one fractal dimension, 
then to remove features of less than two fractal 
dimensions, and finally to remove features of less than 
three fractal dimensions from the original three- 
dimensional data set as the morphological skeleton is being 
formed. Thus, the cascade of filters is used directly with 
a tomographic image to isolate anatomical structures from 
surrounding tissues to facilitate analysis and review 
thereof . 

By utilizing the residuals of morphological erosion 
and opening, the morphological skeleton is formed. This 
process is ideal for processing data with fractal 
dimensional components. For example, the recursive 

formation of the morphological skeleton utilizing 
alternating opening and erosion transforms a 3.4 
dimensional form into .4 dimensional data when a spherical 
structuring element is utilized. 
5 Once the morphological skeleton has been formed via 

recursive development utilizing alternating opening and 
erosion processes, then fuzzy connectivity is utilized in 
the reconstruction of those anatomical structures of 
interest. Reconstruction of anatomical structures without 
0 utilizing fuzzy connectivity results in the loss of 
significant features such as surface textures and 
roughness. These features must be reconstructed from the 
residuals defining the morphological skeleton utilizing 
fuzzy connectivity. The reconstruction of such anatomical 
15 features requires the satisfaction of a fuzzy connectivity 
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criteria such that only those tissue features connected to 
the dimensional features isolated by the sieving process 
are utilized. 

The final result of both the sieving and fuzzy 
connectivity processes is a classification and clear 
visualization of the anatomical structures of interest, 
e.g., tissues and/or tumor pathologies. Additionally, 
quantification of the volume of organs and tumors as well 
as other measurements of interest, such as the diameter of 
arteries and veins, are easily facilitated as a direct 
result of the use of dimensional sieving and fuzzy 
connectivity . 

Connectivity is a mathematical concept which states 
that a set of points is connected if and only if every pair 
15 of points in the set can be connected by a line which is 
contained within the set. The algorithm described in this 
invention generalizes this concept of connectivity to the 
discrete topological grids utilized by a computer to store 
the digital image data by utilizing fuzzy set operators. 
A fuzzy set is itself a generalization of a discrete set by 
defining a function over a set representing degrees of 
membership such that membership varies from zero which 
indicates no membership to one which indicates complete 
membership • 

25 To define connectivity, this algorithm utilizes a 

fuzzy generalization of mathematically defined distances 
between sets as a connectivity criterion. This criterion 
establishes that if two points or two sets of points are 
within a specified distance of one another, then they have 

30 membership to the same set of points. 

The prior art attempted to isolate anatomical features 
from one another based solely upon the intensity of pixels 
within the. three-dimensional data set. The present 
invention facilitates the distinguishing or isolation of 

35 anatomical features based upon such criteria such as size. 



20 
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shape, and intensity of the anatomical feature. Thus, more 
flexibility in designating those features to isolate is 
provided and improved accuracy of such isolation is 
attained . 

These, as well as other advantages of the present 
invention will be more apparent from the following 
description and drawings. It is understood that changes in 
the specific structure shown and described may be made 
within the scope of the claims without departing from the 
spirit of the invention. 

Brief Description of the Drawings 
Figure 1 is an illustration of the recursive 
alternating opening and erosion processes for two 
dimensions utilized to define the residuals of which the 
morphological skeleton is constructed; 

Figure 2 shows the two-dimensional structuring element 
utilized in the process for forming the morphological 
skeleton shown in Figure 1; 

Figure 3 is a chart giving the results of utilizing 
structuring elements of different forms or dimensionalities 
upon images of different forms or dimensionalities; 

Figure 4 shows a representative two-dimensional 
structuring element utilized in the fuzzy connectivity 
restructuring process wherein 2r is the major diameter 
thereof ; 

Figure 5 shows the use of the structuring element of 
Figure 4 to determine that two points belong to the same 
set, i.e., a set of data points defining a desired 
anatomical structure for reconstruction, the two points 
belong to the same set since when one of the points is 
located at the center of the structuring element, the other 
point falls within the bounds defined by the structuring 
element, wherein the dimension d defines the dimension 
between adjacent points such that the points fall within 
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-the set; 

Figure 6 shows the use of the restructuring element of 
Figure 4 to iteratively determine that the set points 
illustrated are contained within a common set; 
5 Figure 7 shows the set support function which defines 

the degree of fuzzy membership for a given pair of points, 
which is determined by the modified Hausdorff metric for 
those points. 

Figure 8 is a block diagram of the conventional 
10 morphological data decomposition and reconstruction 
processes; 

Figure 9 is a block diagram of the morphological data 
skeletonization process of the present invention; 

Figure 10 is a block diagram of the morphological data 
15 decomposition and selective reconstruction processes of the 
present invention; 

Figure 11 is a block diagram of the morphological data 
dimensional sieving decomposition and selective 
reconstruction processes utilizing a three-dimension 
20 exaonple; 

Figure 12 is a block diagram of the morphological data 
reconstruction from skeleton process of the present 
invention ; 

Figure 13 is a block diagram of the morphological data 
25 decomposition and selective reconstruction process of the 
present invention; and 

Figure 14 is a block diagram of the fuzzy logic 
process of the present invention. 

30 Detailed Description of th e Preferred Embodiment 

The detailed description set forth below in connection 
with the appended drawings is intended as a description of 
the presently preferred embodiment of the invention, and is 
not intended to represent the only forms in which the 

35 present invention may be constructed or utilized. The 
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description sets forth the functions and the sequence of 
steps for constructing and operating the invention in 
connection with the illustrated embodiment. It is to be 
understood, however, that the same or equivalent functions 
5 and sequences may be accomplished by different embodiments 
that are also intended to be encompassed within the spirit 
and scope of the invention. 

The data dimensional sieving and connectivity 
methodology of the present invention is illustrated in 

10 Figures 1-14, which depict a presently preferred embodiment 
of the invention. 

Referring now to Figure 1, the recursive development 
of a morphological skeleton utilizing alternating opening 
and erosion process is shown utilizing a two-dimensional 

15 geometric construction, i.e., a square, for purposes of 
illustration. Although a two-dimensional example is 
provided herein, for purposes as illustration, those 
skilled in the art will appreciate that use of the present 
invention in medical imaging typically requires the 

20 recursive use of a three-dimensional structuring element 
such as a sphere, a two-dimensional structuring element 
such as a surface, a one-dimensional structuring element 
such as a curve, and a zero-dimensional structuring 
element , i.e., a point • 

25 After the first opening process, a square 101 having 

the corners removed therefrom is defined. An octagon 100, 
as shown in Figure 2, is utilized as the structural element 
for this example. The corners 102 are the residuals of the 
opening process for the original square. Each time an 

30 additional erosion and opening process is performed, 
progressively smaller squares 102, 103, and 104 is formed. 
After each recursive erosion and opening process, 
additional residuals 102 are defined. After the last 
erosion process is performed, the square is completely 

35 eliminated and the collection of residuals defines the 
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desired morphological skeleton 106. 

Dilation and erosion are defined as follows: 
erosion, 

(gef ) (x)=rain{g(x+d)-(f.(d)-f JO) ) >, 

* d in E 



dilation, 

(gef J (x)=max{g(x+d) + (f J-d)-f,(0) )} 

* d in E 

= - (- g e fj(r) 

Structuring element, f , and image function, g, defined over 
domain of definition for f, E, 

f.(x) = a f(i) ^ 

d=inf{ d inf (Xag)=sup(X„g) ,oe3l+}, 
ag = g(x/a),3l+ denotes the real numbers iO 
Se„(g,X) = {Xeag)-(Xeag)+eg, +0-limeioe 

Alternately for a black skeleton the extensive 
operations of dilation and closing are performed. 

b = inf{ ol inf (X«g)=sup(X«g) ,o6a+}, ag - g(x/o), 

Se_„(g,X) = (X®og)-(Xeag)+e9, a^O 

For digital raster formats of pixels or voxels, d is 
limited to the integer domain Z of the data and e is equal 



to 1. 

I 
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defined as a single dilation step followed by a single 
erosion step. 

By decreasing the size of the structuring element 100, 
smaller residuals 102 are obtained and the resolution of 
the morphological skeleton is increased. 

This morphological skeleton contains all of the 
information contained in the original image. The original 
image can be reconstructed from the morphological skeleton 
by reversing the recursive development process, i.e., by 
substituting dilation and closing for erosion and opening, 
respectively. Thus, by performing a series of dilations 
and closings, instead of the openings and erosions 
performed previously, the original three-dimensional data 
set is obtained from the morphological skeleton. 

In forming the morphological skeleton 102, data 
dimensional sieving is performed such that anatomical 
structures having various dimensionalities are separated 
from one another in a manner which isolates them and makes 
them identifiable via computational methodology. Thus, 
according to the methodology of the present invention, 
those anatomical structures having a fractal dimensionality 
of less than one dimension are separated from those 
anatomical structures having a fractal dimensionality of 
less than two dimensions, both of which are separated from 
anatomical structures having a fractal dimensionality of 
less than three dimensions. 

A desired anatomical structure which has been so 
isolated and identified can then be reconstructed by 
reversing the recursive morphological skeleton development 
sequence described above utilizing only the data points 
associated with the selected anatomical structure. 
However, merely reconstructing the desired anatomical 
structure results in the loss of significant features such 
as surface textures and roughness. 

Thus, in order to preserve such significant features. 
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it is necessary to utilize fuzzy connectivity during the 
reconstruction process. The use of fuzzy connectivity 
assures that all of the data points associated with the 
anatomical structure are utilized in the reconstruction 
process. According to the present invention, fuzzy 
connectivity defines the entire data set for the desired 
anatomical structure by utilizing a modified Hausdorff 
metric, wherein connectivity is defined by the size and 
shape of the structuring element. 

For example, the structuring element is first centered 
upon a seed pixel by the operator. The seed pixel is one 
which the operator knows is a part of the anatomical 
structure for which reconstruction is desired. All other 
pixels contained within the volume defined by the 
15 structuring element are then considered to be a part of the 
anatomical structure being reconstructed. This process is 
then repeated for each new pixel within the data set until 
no additional new pixels are found. Although, as in the 
formation of the morphological skeleton, many different 
20 sizes and shapes of structuring elements are suitable, 
those generally spherical in configuration are preferred. 

A series of different structuring elements may be 
utilized in either of the formation of the morphological 
skeleton or the reconstruction process, as desired, so as 
25 to achieve a desired effect. 

As mentioned above, connectivity is a mathematical 
concept which states that a set of points is connected if 
and only if every pair of points in the set can be 
connected by a line contained in the set. The algorithm 
30 described in this invention generalizes this concept of 
connectivity to the discrete topological grids of computers 
and digital image data with fuzzy set operators. A fuzzy 
set is itself a generalization of a discrete set by 
defining a function over a set representing degrees of 
35 membership from no membership as represented by a zero to 
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complete membership as represented by a one. This 
algorithm utilizes convex fuzzy membership, as shown in 
Figure 7, functions defined over convex set supports. 

To define connectivity, this algorithm uses a fuzzy 
generalization of mathematically defined distances between 
sets as a connectivity criterion. This criterion 

establishes that if two points or two sets of points are 
within a specified distance of one another, then they have 
membership to the same set of points. To more precisely 
define this concept of connectivity, the neighborhood of 
points and the data must be defined. 

As shown in Figure 7 , convexity implies that a line 
fixed between any two points on the curve of the function 
must lie on or below the graph of the function: 

k f(a) + (1 - A) f(b) ^ f(Aa + (1 - A)b), 

erosion, 

(gef J {x)=min{g(x+d)-(f^(d)-f^(0) )}, 

d in E 

dilation, 

(gef J (x)=max{g(x+d) + (f^(-d)-f^(0) )} 

* d in E 

= - (- g e fj(r) 

Structuring element, f , and image function, g, defined over 
domain of definition for f, E, 

f,(x) = a f(§) ^ 

Minimum Function Maximum function 

(g A h) (x)=min{g(x) , h(x)>. (g V h(x)=max{g(x) , h(x) } . 

Based on the previous definitions a measure of distance 
between sets or points g, h can be defined. This metric is 
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used to define points or sets within this distance to be 
fuzzy connected. 

2.4.1 Definition: Modified Hausdorff function 

metric 

df(g,) = minj o such that max ( (gAh) (x) ) 

< min( ( (g®f.) A hef.) ) (x) ) ) . 

X in E * ' 



Referring now to Figure 3, a chart showing the result 
of utilizing a structuring element of a particular form or 
dimensionality on an image of a particular form or 
dimensionality is shown. The chart includes structuring 
15 elements of point, segment, disk, and sphere form and 
images of point, curve, circles, and volume form. As shown 
in the chart, utilizing a structuring element defined by a 
point, for example, in the processing of a curve according 
to the methodology of the present invention, yields a 
curve. Similarly, utilizing a segment in the processing of 
a curve yields a curve and utilizing a disk or sphere in 
the processing of a curve provides a null product, since a 
two-dimensional disk or a three-dimensional sphere cannot 
be utilized to process a one-dimensional curve. 

Referring now to Figures 4-6, the use of a two- 
dimensional example of a structuring element and the fuzzy 
connectivity reconstruction of a desired anatomical 
structure is shown. With particular reference to Figure 4, 
the structuring element 200 shown comprises an ellipse 
having a major diameter of 2r. Those skilled in the art 
will appreciate that various other shapes are likewise 
suitable for use as a structuring element. 

Referring now to Figure 5, use of the structuring 
element to determine if two points are within a common set 
35 is shown. This is accomplished by placing the structuring 
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element 202 around one of the points 210 of interest and 
then determining whether or not the second point of 
interest 212 lies within the boundary of the structuring 
element 202. As shown, the second point 212 does lie 
5 within the boundary of the first structuring element 202. 
In order to find additional points which are part of the 
common set of points, which define the anatomical structure 
of interest, this process is repeated by placing a 
structuring element 204 around the second point 212 in 
10 order to determine if any points lie within the boundary 
thereof . 

With particular reference to Figure 6, this process is 
repeated so as to define all of the points which belong to 
a common set of data points which define the anatomical 

15 structure of interest. Structuring element 202 formed 
about point 210 defines point 212 as being included within 
the data set, structuring element 204 formed about point 
212 similarly defines point 210 as belonging to the common 
data set, while structuring element 206 formed about point 

20 212 defines point 214 as belonging to the common data set. 
Thus, all points which lie within the boundary of any 
structuring element at which a point within the data set is 
formed at the center thereof, also are members of the 

common data set. 
25 Each point so defined to be within the data set is 

assigned a fuzzy membership number between zero and one, 
depending upon the distance between adjacent points, as 

discussed above. 

Thus by utilizing fuzzy connectivity, the set of all 

30 data points defining a particular anatomical structure of 
interest are defined such that surface details of the 
anatomical structure, such as surface smoothness thereof, 
are maintained during the reconstruction process and are 
thus included in the reconstructive anatomical structure. 

35 Referring now to Figure 8 , an overview of the standard 
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morphological decomposition and reconstruction process is 
shown. According to contemporary methodology, an input 
data array 300 is skeletonized 302 so as to form skeleton 
304. Skeleton 304 is then reconstructed 306 so as to 
5 provide the original image 308. This process is used in 
various different data analysis, compression, and data 
signal processing applications. 

Referring now to Figure 9, morphological data 
skeletonization according to the present invention is 
10 shown. Morphological data skeletonization is a recursive 
process wherein erode image n 320 subjected to erosion 
322. The product of erosion is then subjected to dilation 
323 and in parallel subjected to erosion 324. The product 
of erosion 324 is erode image n+1 326 which then becomes 
15 new erode image n 320 and is iteratively processed. The 
product of dilation 323 is subjected to subtraction 325 
with respect to erode image in 320 so as to form skeleton 
327 which is then subjected to addition with full skeleton 
304. 

20 Referring now to Figure 10, morphological data 

decomposition and selective reconstruction according to the 
present invention is shown. Input data array 300 is 
subjected to skeletonization so as to form skeleton 304. 
Skeleton 304 is used for the selection of a region of 
25 interest 310 so as to form edited skeleton 312. Fuzzy 
connectivity 314 is applied to the edited skeleton 312 to 
form the edited image 316. 

Referring now to Figure 11, a three-dimensional 
example of the process of morphological data dimensional 
30 sieving, decomposition, and selective reconstruction is 
shown. Input data array 300 is skeletonized 342 wherein 
a three-dimensional kernel or structuring element 
configured as a sphere, for example, is utilized in the 
skeletonization process. The skeletonization 342 results 
35 in the formation of a skeleton 343 having less than three- 
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dimensional features. This skeleton is then subjected to 
skeletonization 34 4 utilizing a two-dimensional kernel or 
structuring element configured as a facette. This two- 
dimensional skeletonization process 344 results in a 
5 skeleton having less than two-dimensional features 345 • 
This skeleton having less than two-dimensional features 34 5 
is then subjected to skeletonization utilizing a one- 
dimensional kernel or structuring element 346 so as to 
provide a skeleton having less than one-dimensional 

10 features 304 . 

Referring now to Figure 12, the process of 
morphological data reconstruction from a skeleton without 
the use of fuzzy connectivity is shown. As discussed 
above, such reconstruction results in the loss of 
15 substantial surface detail. Using reconstruction n 

350, dilation 352 is performed so as to produce dilate 
image n 353, dilate image n 353 and skeleton n 354 
are added 356 and the process is iterated by providing the 
added images as recon n 350. 
20 Referring now to Figure 13, the process of 

morphological data decomposition and selective 
reconstruction of the present invention is shown. Recon n 
360 is subjected to dilation 364 so as to produce dilate 
image n 368 and seed image n 366. Seed image n is 

25 subjected to fuzzy connectivity criteria 370 with skeleton 
n 362 so as to produce edited skeleton 378. Dilate image 
n 368 is combined 380 with edited skeleton 378 so as to 
produce a new recon n 360 and the process is iterated. 
Referring now to Figure 14, the use of fuzzy 
30 connectivity according to the present invention is shown. 
A seed image n pixel 400 and the skeleton n 402 are 

operated upon by fuzzy logic 404 utilizing pixel fuzzy 
logic measure 406, i.e., the selective structuring element, 
so as to provide pixel fuzzy measure update 408 and set 
35 fuzzy connected pixels 410. The use of fuzzy logic in 
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this manner is described in detail in "Analysis and 
Segmentation of Higher Dimensional Data Sets With Fuzzy 
Operators for Representation and Visualization" and 
published in w^uro and F n^-^iv Systems; — Emergent Science of 
Tntelliaent Computing , by Mitra, Gupta, and Kraske, 
published by SPIE Press, 1994, ISBN 0-8194-1566-9, provided 
herewith and forming a part of this patent application, the 
entire contents of which are hereby incorporated by 
reference. 

Provided below is a list of the symbols utilized in 
the math equations in this patent application: 
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— Module space of integers of dimension n. 
Z+ — Half space of integers greater than zero. 

— Euclidean space of real numbers of dimension n. 

E — Space domain of which the image is defined over, a Borel field, locally compact, 

Hausdorff and separable, t>pically a compact or SR^ space. 
CI — Set of aU directionf: from the center of the unit sphere angles of 4tz vieradlans. 
CO — Arbitrary angle within the range of the unit sphere 47c steradians. 
^ — A point in E of coordinates (x,y,z). 

r — A limit point in E of coordinates (x,y,z). 

OC, d, b — Points on the real line, 91, or integer domain, Z. 
d> b — Limit point on the real line, 91, or integer domain, Z, 
d(rl j2) — Euclidean distance between points rl and r2. Vcx.l-x2)2+(yi-y2)-+(zl-z2)- 
B — A stiucniring element set, compact, convex, symetrical to inversion, 

with unit radius, .5 • jB] = 1. 
|B| — Diameter of a set : sup { d(rl,i2) 1 rl,r2 e B } 

Bf — A structuring element set with center translated to point r in the space E. 

CCB — Scale of the set B by factor a uniformly in all directions relative to the origin 

of the set coordinate system of definition. 
X — A set in the space E or an image defined over E. 

— complement of the set X in the space E. 
0 — The null or empty set. 

© — Image algebraic erosion, Minkowski set subtraction. 

ffi — Image algebraic dilation, ^*linkowski set addiUou. 

X — The closing of a set in the space E or an image defined over E with the 

structuring element B. (X © B) © B. 

Xg — The opening of a set in the space E or an image defined over E with the 

structuring clement B. (X 0 B) © B. 
V(0 _ The measure defined over the space E of functional distributions mapped over 
the sets of E. A Lebesque measure over continuous space or the sum of voxel 
intensities over the integer lattice of digital imagery. 

A — The Infimum of the functional mappings defined over the cuclidean space, 
for integer lattice of digital imagery this corresponds to the minimum of the 
functional intensity mappings of each voxel. 

V — The supremum of the functional mappings defined over tiie euclidean space 
for integer lattice of digital imagery this corresponds to the maximum of the 
functional intensity mappings of each voxel. 

/ — The set difference of binary sets in the euclidean domain or functional 

difference of mappings over the euclidean domain. 
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Provided below is a program listing utilized in the 
practice of the present invention upon a Pixor computer: 
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chmap -i 



VolClaGGify -pat LOWERY -map M mattetable R loweryvent -roi 40 27q ocn ^ \ 
-in 60 -first 4 -name /usr/wolf /data/loweryvent ^ -A Z \ 

ffchmap -X 

gVolclaccify -pat LOWERY -map M mattetable R loweryvent -roi 40 279 0R^ > ^ 
g 7in 60 -first 4 -name /usr/wolf /data/loweryvent ^ ^ 

#chmap -1 

# VolClaseify -pat LOWERY -map M mattetable R lowerybrain -roi 4 0 279 -^a t 
^ -A Z -in 60 -first 4 -name /usr/wolf /data/lowerybJain ^ 38 261 -1 
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#include "DBM sizes. h" 

^include "VOX^iitnagelist - h" 

^include <merge . h> 

^include <math . h> 

^include <vv__impl . h> 

# include <pixar /video , h> 
ttinclude <pixar/chap/pw. h> 

# include <videor eg . n> 
#include <varargs . h> 
#include <sys/time . h> 
#include <sys/etat .h> 
# include "VOXAR . h" 

#define PIRL VOL NULL 39 
Sdefine PIRL"WI>rNULL 40 

#define pirl3tbl3;;null 41 

#define INV DBL2PXL(1.2) 

^define mod{a, b) (((a) < 0) ? {(b) + (a) % (b) ) : ((a) % (b) ) ) 

extern RGBAPixelType RGBAzero, RGBApoint, RGBAquart, RGBAhalf, 

RGEAone, RGBAmsLX, RGBAinv, RGB An 

extern COLOR_MATRIX_type CHANzero, CHANident; 

extern VIDEO *ThePirlVideo ; 

extern double *make_laplace () , sqrt () ; 

PirlError voxVol Skeleton (pv, vol_tbl , thintable , thinscale , thincnt) 

PirlPV ♦pv; 

vol table type *vol tbl ; 

RGBSPixelType * thintable ; 

Rouble yg^^g,^gj^^^^^®^^^^f^ij^„jLdth. thinheight, thindepth, pvdepth, p, threshold; 

static^^RGB^ixelTyp"' thintablel [4096] * arowtable r409fi3 , mawVi-ablo (aos) 
a>uciLit; ^Xx-iFV "sfcSl, *ait£, erod, smoothpv; 
static double *k, *ka; 
static TT-5/ sum* value; 

static char FuncName[) = »'VOXVolSkeleton» ; 
FILE *fp; 

if ((thincnt & 0x3F) 1= thincnt) thincnt « 0; 
fprintf (stderr, "function Vs\n" , FuncName) ; 

if(lT>v 11 •*Dv II ivol tbl I Ivol tbl->vol_link J) tthxntable) 
{ fprintf (stderrT "Vs cannot"load a null volume", PMncName) ; 

CHECK (PIRL_VOL_NULL) ; 

fvdepth «» CountSlices{vol_tbl) ; 
printf (stderr, "%s thin^ecale - %f" FuncName, thxnscale) ; 
'*/ if (thinscale 0. | thinscale > 1,5) thinscale « 1.5; 
if (thinscale > !•) TT=3; else TT-5; . 

if(5(k =. make 3D sep(k, TT, vol tbl, thinscale))) 

fprintf (etderr, "kernel cannot be defined"); 
dump_Bep_kernel(k, TT, stderr) ; 

fprintf (stderr, **%s volume x - Vd y - Vd z • %d\n", FuncName, 
(*pv) ->xsize, (*pv) ->ysize, pvdepth) ; 

thinwidth ^ (int) (thinscale * (double) ( {*pv) ->xsi2e) ) ; 
thinhe.lghh • (int) (thinscale *• (double) { (*pv) -^ysize) ) ; 
zspace = (double) (vol_tbl->sl ice. pix) ; 

/* th?ndtpth i^uSt) '(thinscale * {*^*=»^^?-^l^rP'^ir??iv''U4 1 «o« , 

♦/ zspace /« vol tbl->3lice_pos . z - vol_tbl->vol_link->slice_pos.z; 

zspace (zspace < 0.) ? -zspace : zspace; 
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thindepth = (zspace > 1.) 

? pvnp.pth : (int) (sspacc* (aouble) (pvdepth) ) ; 
CHECK ( 

l^ll^^i^^^^ (*pv}.>xsize, (*pv)->y3ize, pvdepth)); 

PirlCreatePV(&erod, ( *pv) - >xsi2e, (*pv) --.ycize, pvdepth) ) ; 

fprintf (stderr, "%s volume x = %-d y = %d z = %d\n" FuncNam^ 
smoothpv->xsize. smoothpv* >ysize , 6moothpv->zsize) ) 

fprintf (stderr, "%s thin width - hei^^^ T '^^^^P^^* - ^^<i\n% FuncName, 

— thxnwidth, thinhexght, thindepth); 

for(p = 0; p < 2048; p++) 

if (thintable [p] .Alpha) thintabiel tp) = RGBAcne; 

thintabiel fp] « RGBAzero 

)ca-= (double *) malloc {27*sizeof (double) ) ; 
p-27; while(p--) ■ 
SUm»0; g-27; 
threshold = DBL2 
if (threshold <i- 8*DHL2PxL'(*ka) ) 

fprintf (stderr, ">\c!\t r.h-reshold is too LARGn\n"); 
rprxntt (stderr, "threshold minimum ^ Vd\t", 8*DBL2PXL ( *ka) i • 
fprintf (stderr, "kernel sum = %d\t", sum); ^^^^-^^^ ^a; ; , 

fp=f open ( "kernel . file" , "w") ; 
dump_kernel (ka, fp) ; fclose(fo) ; 
fprxntf (stderr , "threshold = %d\h", threshold); 
for{p = 0; p < sum-threshold; p++) 
for(p = sum- threshold; p < 2048; p++) 



"i-n / Jca(p]= k(p/9) * k[(p % 9)/3] * k[p^3] 
>p5akl^I]1-l; ^""^ DBL2PXL(ka[p]); 



0, 0) ) ; 



for(p = 0; p < threshold; p++) 
for(p ^ threshold; p < 2048; p^-*) 

ShowTree ( PV7StackRoot , FuncName) ; 

CHECK ( VOXFas tDumoVolume (pv , NULL) ) ; 
CHECK ( VOXThreshVblume (pv, thintabiel , 
if(thinscale > 1.) { 
CHECK ( VOXConvVolKernel (pv, pv, ka) ) ; 
CHECK ( VOXThreshVolume (dv, growtable, 0, 0)); 

CHECK ( VOXCopyVolume (pv] fcerod) ) ; 
p = thincnt; 

while (thintable fitfic p-*) { 

f print f (stderr, "%s thin count = %d\n", FuncName, p) 
CHECK ( VOXCopyVolume ( ficerod, tamoothpv) ) ; 
CHECK ( voxconWolKernel (tsmoothov, &ercd, ka) ) ; 
CHECK( VOXThreshVolume (fiierod, masktable, 0, 0)); 
CHECK ( VOXConvVolKernel (fiterod, 6cSmGOthpv, ka) ) ; 
CHECK { VOXThreshVolume ( fcsmoothpv , growtahl ft , 0 , 0 ) ) 
CHECK ( voxsubvolume (&smoothpv, «terod) ) ; 
CHECK ( VOXSubVolume (pv, &smoothpv) ) ; 

. ) 

CHECK ( VOXAddVolume ( &erod , &smoothpv) ) ; 
p - thincnt; 

while (thintable p--) { 

fprintf (Btderr, "%s thin count « %d\n", FuncName, p) 
CHECK( VOXConvVolKernel (ficerod, &smoothpv, ka) ) ; 
CHECK ( VOXThreshVolume (ismoothov, growtable, C, 0) ) 
CHECK ( VOXSubVolume (tsmoothnv, £cerod) ) ; 
CHECK ( VOXMul t Volume (fitsmoothov, pv) ) ; 
CHECK ( VOXAddVolume {&erod, ^smoothov) ) ; 
CHECK ( VOXConvVolKernel (&erod, terbd, ka) ) ; 
CHECK( VOXThreshVolume (&erod, growtable, 0, 0)) 



1^ ^."i^sktable tp] = RGBAzero; 
masktabletp) - RGBAone; 

growtable tp) RGBAzero; 
growtable tp) = RGBAone; 



p » thincnt ; 

if (thinscale > 1.) 

while (thintable 



) ( 
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forintf (etderr , "ts thin count - %d\n" , FuncName, p) / 
CHECK ( VOXConvVolKernel (Scerod, &smoothpv, ka) ) ; 
CHECK ( VOXThreshVolume ( &Bmoothpv, growtable , 0,0)); 
CHECK ( VOXSubVolume(&:smoothpv, &erod) ) ; 
CHECK ( VOXMultVolume ( &smoothpv, pv) ) ; 
CHECK { VOXAddVolume (&erod/ &smoothpv) ) ; 
CHECK ( VOXConvVolKernel (terod, &erod, ka) ) ; 
CHECK( VOXThreshVolume (&erpd, growtable, 0, 0) ) ; 

CHECK { VOXFastReadVolume (pv, NULL) ) ; 

CHECK { VOXMultVolume (pv, terod) ) ; 

CHECK ( PirlDeletePW(&:smoothpv) ) ; 
CHECK ( PirlDeletePW(ficerod) ) ; 

return PIRL NO_ERROR; 
error : reuuiii PirlLastiirr ; 

} 

PirlBrror VOXThinVolume (pv, vol_tbl, thintable, thinscale, thincnt) 

PirlPV *pv; 
vol table type *vol_tbl; 
RGBSPixelType * thintable ; 
double thinscale; 

{ register thinwidth, thinhexght, thindepth, pvdepth, p, thresholds- 

register double zspace; , r , . , , , 

static RGBAPixelTyoe thintablel [4096] , growtable [4096] , masktable [409 
static PirlPV origpv, newpv, pvresize; 
static double *k; 

static TT=5, vent flag=0; , . , 

static char FuncName [] « "VOXThinVolume" ; 

if (thincnt < 0) { thincnt- -thincnt; vent_flag«0; } 

else vent__f lag-0 ; 

if ((thincnt t Ox3F) != thincnt) thincnt «= 0; 

fprintf(stderr. "function %s\n% FuncName); ^. ^ 

ifdpv II !*pv II !vol tbl II !vol_tbl->vol_link Jl ithintable) 
7 f nr-inhf (stderxT *'kc cannot load a null volume", FuncName); 

CHECK ( PIRL_VOL_NULL ) ; 

if (thinscale > 1.) TT=:3 ; else TT«5; 

pvdepth =: CountSlices (vol tbl) ; 

lf(pvdepth > (*pv) ->2sizeT PY4^P^Ji*"'„ 
/* fprlntf i[stderr, »%s thin scale = %f% FuncName, thinscale) ; 

*/ if (thinscale <= 0, \\ thinscale > 2.) thinscale =2,; 

if{!(k ^ make 3D 6ep(k, TT, vol_tbl, thinscale))) 

fprintf (stderr, "kernel cannot be defined"); 

dump_sep_kernel(k, TT, stderr) ; 

fprintf (stderr, "%s volume x = %d y = %d z = %d\n" , FuncName, 
(*pv) ->xsize, (*pv) ->ysize, pvrippth) ; 

thinwidth - (int) (thinscale * (double) ( (*pv) ->xsi2e) ) ; 
thinheight = (int) (thinscale * (double) (*pv) ->ysi2e ) ; 
/* thindepth = (int) (thinscale * (double) ( (*pv) ->2si2e) ) ; 

zspace = (double) (vol_tbl->slice .pix) ; 

zspace /- (1000.); ^ .... . 

*/ zspace /= vol tbl->slicejpos . z - vol_tbl ->vox_link- >slice_pos . z ; 

zspace = (zspace < o.) ? -zspace : zspace; 
thindepth « (zspace >1.) • ^v.x\ 

? pvdepth i (int) (zspace* (double) (pvdepth) ) ; 

/ * CHECK ( 

PirlCreatePV(torigpv, (*pv) ->x3ize, (*pv) ->ysize, pvdepth) ) ; 
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* / CHECK ( 

rix-lci:caLc:FV(dcnewpv, I'pv) ->x5i2e, (*pv) ->ysi2e, pvdepth) ) ; 

fprintf (etderr, "Vs new pv width « %d height - kd depth := %d\n", FuncNam 

fprintf (stderr, "Vs thin width » %d height'T''id''dijth S^vI\;;"r-uncNam^ 
^/ thinwidch, thinheight, thindepth)"; ' 

forCp = 2048; p < 4096; p++) • thintablel (p] . RGBAzero; 

for{p 0; p < 204B; p++) 

if (thintable [p] .Alpha) thintablel [p] « RGBAone; 

else thintablel [p] = RGBAzero 

threshold « DBL2PXL ()c [11 ♦ *k * *)c) ; 
threshold - 256; 

fprintf (stderr, "threshold - Vd\n", threshold); 

for(p ^ 0; p < (2048-threshold) ; p-i.+ ) inas)ctable tp] RGBAzero; 

for(p =t (2048-threshold); p < 2048; p++) masktabla [p] - RGBAone; 

for(p = 0; p < threshold-20; p++) growtable [p] <= RGBAzero; 

tor{p = threshold; p < 2048; p-^+) growtable tp] RGBAone; 

CHECK ( VOXCopyVolume (pv, &newpv) ) ; 
ShowTxeti (PWScackKoot , FuncName) ; 

CHECK( VOXThreshVolume (tnewpv, thintablel, 0, o)); 
p - thincnt; 

while (thintablel p--) { 

fprintf (stderr, "%s thin count = %d\n", FuncName, p) ; 
CHECK ( VOXSepConvVol ( fcnewpv , &nswpv , k , TT) > ; 
CHECK( VOXThreshVolume (&newpv, masktable, 0, 0) ) ; 

if (vent flag) p « thincnt/2; 

else p = thincnt; 

while (thintablel && p--) { 

fprintf (stderr , "%s thin count = Vd\n", FuncName, p) ; 

CHECK ( VOXSepConvVol (&newpv, &newpv, k, TT) ) ; 

CHECK ( VOXThreshVolume ( toewpv , growtable ,0,0)); ^ 

CHECK ( VOXMult Volume (pv, tnewpv) ) ; 

CHECK ( PirlDeletePW (&newpv) ) ; 

return PIRL NO ERROR; 
error: return PirlEascErr; 
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#include 

include 
4finclude 

include 
# include 
S include 
(J include 
# include 
# include 
^include 
# include 
^include 

ttdefine INV 



"DBM sizes. h" 
''VOX^Iimageliet.h" 
-tfmerge .h> 
<macn.h> 
<w_impl ,h> 

<pixar/video . h> 

<pixar/chap/pw.h> 

<videoreg.h> 

<varargs.h> 

oys/txme.h> 

<sys/stat .h> 

"VOXAR.h" 



DBL2PXL(1,2) 



#define PIRL VOL NULL 
#define PIRL~WIN"KXJLL 
#define PIRL_TBL~NULL 

#define mod{a, b) 

extern RGBAPixelType 



39 
40 
41 



(((a) < 0) 



? ((b) + (a) % (b)) 



((a) V (b))) 



RGBAzero , RGBApoint , RGBAquar t , RGBAhalf , 

RGBAone, RGBAmax, RGBAinv, RGBAn 



extern C0LOR_MATRIX_type 



CHANzero , CHANident ; 



extern VIDEO 
extern double 
PirlError 



PirlPV 

vol table_type 
douole 

RGBAPixelType 
char 



*ThePirlVideo; 

*make_laplace ( ) , sqrt ( ) ; 

VOXProiect (in_pv, vol__tbl, perm, anqle, range, incr, 
gradtable , reconcable , f _name , beam^vol , beamtable ) 

*in_pv, *beam vol; 
♦voTTtbl ; " 
angle, range, incr; . 

*gradtable , *recontable , *beamcable ; 
*f_name; 

PirlPW srcpw, ramp_pw, grad__t, grad_m, grad_b, beampw, 

destpw, pt pw, scratchpw, prevpw, gradow, reconpw, 
shade , shaclesub , buf j)w , temp^pw ; 
RGBAPixelType *bgtable, *tumortable, 

RGBAchanl= { 0x300 , 0 , 0 
double *k, factor; ^ 
coord3 D_dbl g , resize; 

char *line_8tr, FuncName t] « "VOXProject" , 

filename [12 8] ; 
rec_width, rec_height, 

cnt , num^^cut , 1 ; 
shadeh , shadew, pw^height , pw^width, xof f , yof f ; 
fd; ^ 

struct { short multichannel , width, height , minval , max 

} FileHeader; 

line str « malloc(512); 

bgtaEle* (RGBAPixelType *)malloc (2049*aizeof (RGBAPixelTyoe) ) ; 
if (beamtable (800] .Red) 

?et_mappinQ table (bgtable, "beamgrad") ; else bgtable-tNULL; 
printr (stderr, "executing %s\n", FuncName) ; 
iidin pv II !*ih pv || !vol tbl !vol tbl->vol link) 
T "^CHEtK( PTRL_V6l_NULE) ; ' * T 

if ( ! recontable I | Ibeamtable) 

{ fprintf (stderr , "%s input table is null\n", FuncName); 

J goto noerror ; 

if (! gradtable) gradtable ^ recontable; 
range angle; 

fprintf (stderr, "%s angle : initial- %e final^^ %e increment* %e\n", 

FuncName , angle , range , 



static 

static 

static 
static 
static 

static 

static 
int 

static 
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factor » (double) (vol tbl->slice .oix) /lOOO . • 

factor /= (vol tbl->sTice__pos,2 vol^tbl- >vol link->slice x>o3 2) • 
if (factor < O-T factor = -ractor; " ^±x^^_poB .z) , 



>2Sl2e ; 

>2si2e) 



switch (perm & WXYZ) f 
case WX2 : ^ 

num_cut -= cnt « (*in_pv) ->ysize; 
pw_width « (*in pv) ->xsize; pw height « 
rec width - rec^height «r ( ( *in_pv>~>xsize < 
? ("'"in pv) ->zsize : (*in_pv) ->xsize; 
g•^.= 1.; g.y = factor; g.z ^ 1.; 
resize. X = 1.; resize. y = factor; resize. z 

break; 
case WY2 : 
num_cuc = cnt {*in_pv) ->x3i2e; 

pw_width « (*in pv) ->ysi2e; pw height = (*in pv) ->2si2e ; 

width « rec_height = ( (*in_pv) ->ysize < (*in~pv) ->2size) 
? (^in pv) ->2si2e : (*in_pv) ->ysi2e ; 

g.x = 1.; Q.V « factor; O . !?: = T ; 

resize. X = 1.; resize. y « factor; resize. 2 = 1.; 
break; 

default : perm « WXY; 

case WXY : 

num_cut « cnt « (*in_pv) ->2si2e; 

pw_width = (*in_pv) ->xsi2e; pw_height « ( *in ov) - >ysize • 

rec height « rec^width « ( (*in_pv) ->xsize < (*in"pv) ->2size) 
(*in_pv) ->2si2e : (*in_pv) - >xsize; " 



(*in pv) 
(*in2pv) 



1, 



resize. x 
g.x « 1. 



1. ; 
g-y 



g.x 
shadeh 

XOff e 



1. 



resize 
= 1.; g- 

g.y ' 



f acto>-; 
f acuor; 



resize . 2 
} 



1. 



« shadew = (pw width 
(£hadew-pw_wid'Eh) /2 ; 



g.z 1.; 

pv/_height) ? pw_height : pw width; 
yoff = (shaden-pw_heicht)72; 



PirlCreatePW ( &beatnr>w, 
PirlCreatePW(&buf pw, 
PirlCreatePW ( tdes'E^pw, 
Pir iCreat ePW ( &grad_Tn , 
PirlCreatePW ( fitgrad_t , 
PirlCreatePW (fitgradpw, 
PirlCreatePW ( &ramp_jpw. 



pw width, pw height, 0)); 
{5'*"rec widthT/4, (5*rec he 
pw_widch/4, pw^height, Z) ) 
pw_width/4, pw_height, 0) ) 
pw_width/4, pw_height, 0) ) 
pw width, pw_height, 0)); 
snaaew, shaaeh, 0) ) ; 



CHECK ( 

c:heck( 

CHECK ( 

check ( 

CHECK ( 
CHECK ( 
CHECK ( 
CHECK ( 
CHECK 
CHECK ( 
CHECK ( 

PirlCreateSubPW (&shade, tshadesub, xoff, xoff+pw width-l, 

yoff, yoff+pw Height- 1, 
PirlCreatePW (ficprevpw, pw_width, pw_heigKt, 0)) 
PirlCreatePW (ficreconow, rec^width, rec height. 



ght)/4, 0)) 



PirlCreatePW (6ctemp_pw, shadew, shadeh, 0)) 
PirlCreatePW (&scratchpw, pv/_width, pw height, 
PirlCreatePW (icshade, shadew, shadeh, TS)); 



0)) 



CHECK ( 
CHECK ( 



OxF) ) 



fprintf (scderr, 



"Vs windows allocated recon v/idth 
FuncName, rec_width, rec_height) ; 



0)); 
Vd height 



Vd\n" 



CHECK ( PirlBeginLines () } ; 

CHECK( PirlSetLineAt tribunes (line str, 4., PM MERGE, fitRGBAzero, 

PE_FELTTIP7 1 . , PF_SINC, 1 . ) ) ; 



if (! (k 

else 
CHECK ( 
do 

static 
static 



' make_3D_sep (k, 3, vol tbl, 1.))) 
fprintf (stderr, "kerneT cannot be defined"); 
{' cnt=9; while{cnt--) k(cnt) *- 1.5/ 

PlrlRampY (ramp_pw, &RGEAmax, &RG5Azero) ) ; 



coord3D_dbl light_dir = -1.2, .6, -.7 

coord3D_dbl light_dir - .75, .0, 1.5 

static double line [6] 12] ; 

coord3D_dbl resizeT ; 

extern coord3D_dbl *rotate_vec ( ) ; 



}, light. 



Light; 



liaht 
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0 



static PirlPW subrecon; 
static fcnt=0, slice cnt = 0, 

xoffl, yoffl, xoff2, yoff2; 



switch (perm) 
case WYZ: 



light a*-light_dir; 

light_a.y=» 
light a.2«= 
light a«light2:dir; 

light_a.y= 
light_a.2« 
light_a=light_dir; 

rotate vec {flight ,—&light_a, angle, 

fprintT(3tderr, "light direction x = %e y %e z %e", 

light. X, light. y, liaht.z); 
fprintf (stderr, "rotation angle = %e\n", angle); 



case VVXZj 



case WXY: 



^light^dir . z; 
•light_dir .x; 
-light_dir.z; 
■light_dir ,y; 

WXY) ; 



light_a.xs 
brea)c; 

light_a.x= 
brea)c; 



-ligh 
-liah 



2 . ) -1 . ) *pow (sin (tempang) , 2 . ) 




= yoffl+(int) {recnieight*resizeT,y) -1; 



resizeT « resize; 
if (perm & WZ) 
{ double tempang; 

tempang « angle* (asxn (1 •) /SO . ) 
resizeT-x = resizeT. z = 1.; 
resizeT. y = i . + (pow(resize .y , 
resizeT^y « sqrc (resizeT. y) ; ; 

fprintf (scderr. "light direction x « %e y %e z *e", 
resizeT, X/ resizeT. y, resizeT. z) ; 

xof f 1 
yof fl 
xoff2 
yof f2 

PirlClear (scratchpw, &RGBAzero) ; . 
CHECK( Plaice {*in_pv^scratch|jw^^ ->ssize-i, 

CHECK( VTPutPlanc (♦in pv, acratchpw, (*lii_pv) ->zsizs-2, 

" WXY, WSINGLE) ) ; 
CHECK{ VTGetPlane {* in pv, scratchpw, 0/ 
perm|((perm & WZ) ? VWEFLECTX i 0), WSINGLE)); 
CHECK( PirlMulI (scratchpw, JtRGBAmax) ) ; ^aaonx 
CHECK{ VIPMap{scratchpw, scratchpw, WR, gradtable, 2048)) ; 
CHECK( VTCompact (scratchpw, grad_t) ) ; 
CHECK ( 

VTGetPlane (*inpv, scratchpw, 1, ,„,o^»Tr,T v 

perm I ((perm & VV2) ? WREFLECTX : 0 ), WSINGLE) ) ; 

CHECk{ PirlMulI (scratchpw, &RGBAmax) ) ; ^r^A^w 

CHECK( VIPMap (scratchpw, scratchpw, WR, gradtable, 2048)); 

CHECK ( VTCompact (scratchpw, grad_m) ) ; 

cnt = num cut; slice_cnt « 0; 
CHECK ( PirlClear (reconpw, &RGBAzero) ) ; 
while (cnt--) { double centerx, 

CHECK ( PirlClear (shade, &RGBAzero) ) ; 
CHECK( VTGetPlane (*in r>v, scratchpw, 
((3+slice cnt > num cut) ? num cut-1 
perraj ((perm & WZ) ? WREPLECTO : 0) , 
CHECK ( VTGetPlane (♦beam_vol / beampw, 
((3+slice cnt > num cut) ? num cut- 1 
permT((perm t W2) T WREFLECTX : 0), 
CHECK ( PirlMulr(scratchpw, tRGBAmax) ) , 
CHECK ( VIPMap (scratchpw, gradpw, WR, gradtable 

CHECK^^VIPMap (beampw, beampw, WR, bgtable, 2048)); 
else 

CHECK( VIPMap (beamow, beampw/ VVR, beamtable, 2048)); 

CHECK( PirlAdd( gradpw, beampw)); ^^^taw 
CHECK ( PirlMerge (beampw, gradpw, beampw, MergeOpBELOW 



centery; 



: 2+slice cnt) 
WSINGLET) ; 

: 2+slice cnt) 
WSINGLET) ; 



2048) ) ; 

} 
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CHECK ( PirlMulI (beamow, ficRGBAmax) ) ; 
CHHCK ( PirlMulI (beampw, &RGBAmax) ) ; 

CHECK ( VTCompact (gradpw, destpw) ) ; 

CHECK{ VlPMap (scratchpw, scracchpw, WR, recontable. 2048U . 
CHECK( VTGetPlane(*beam_vol, beamow, ^udi^xe, ^o^e;; , 

( (3+slice_cnt > num cut) ? num cut-1 : 2 + slice cnt) 
perm] ((perm & W2) 7 WREFLECTX : 0), WSINGLET) ; 
CHECK( VlPMap(beampw, beampw, WR, beamtable, 2048)); 
CHECK ( PirlMerge (scratchpw, beampw, scratchpw, MergeOpPLUS, 

PirlDivI (scratchpw, &RGBAhalf ) ; 

CHECK ( PirlAdd (scratchpw, beamow) ) ; 

CHECK ( VIPGradient (grad_t , gra3i_m, destnw, gradpw, 

g.x, g.y, g. z) ) ;' 
CHECK( VIPSurface (prevpw, gradpw, shadesub, 1.2, WRGB;^))- 
CHECK( ViPShade (shadesub, gradpw, shadesub, light. x, liaht v i^cht 5 
if?peim--VVY2^ temp^pw, WG, beamtable, 2048) ) ^-gnt.. 

CHECK ( PirlCircularShif t (temp_pw, -15, -3) ) ; 

if (perm«=WXY) ( ^ 

CHECK( PirlCircularShif t (temp_pw, 0, 15)); 

if (beamtable [800] .Red) ^ / 

CHECK ( PirlMul (beampw, temp r>w) ) ; ^ 
CHECK( PirlShurfle (beampw, *^rrr")); 
CHECK ( PirlMulI (beamow, &RGBAmax) ) ; 
CHECK (PirlAdd (shadesub, beampw) ) ; 

CHECK { Pirl Clear (temp pw, &RGBAzero) ) ; 

centerx » (double) (shaaew/2) ; centery = (double) (shadeh/2} ; 
if (light. X > 0.) line[O][0] -0.; 

?lBe line[0][0] -= (float) (shadew- 1) ; 

xr (light. y > 0.) line[0l[l] = 0.; 

else^ ^ line[0]tl) - (float ) (shadeh-1) ; 

ifdignt.x != 0. light. y I- 0.) { double magx, magy; 

magx « (light. x > 0.) ? light. x : - light. x; 
magy « (light. y > 0.) ? light. y : -light. y; 
i£7magx > magy) 
line CO] [1] « (float) 

(centery - (centerx - (double) (linefO] [0]))* light. v / 
else ' ^ 

line [01 to] «: (float) 

(centerx - (centery - (double) (line [0] [1] ) ) * light. x /light.y) ; 

elee i£(light.x == 0.) line[0](0] = (float) centerx; 

else line [01 [1] = (float) cer 

linefll [01 = (float) centerx; linedl fl] = (float) centery ; 

line[lj[0] += shade->xmin; linetOllO] += 6hade->xmin; 



line [1] [13 += shade->ymin; line [0] [1] +« Bhade->ymin; 

lined] [0] - shade ->xmax-l; linefoHO] - shade- >xmax-l ; 

line(l](l] = shade->ymajc; lineLO][l) » shade->ymin; 

fprintf (stderr, •'line limits x :%e y :%e x ;%e y :%e\n", 

line (03 [0], line [0111], line(l]to], lined] [i]); 

CHECK( PirlPolyLine (shade, 2, line)); 

CHECKC PirlRotate (shade, temp pw, shade, angle, 1., 1., 

centerx, centery, 1, 4;;; 
CHECK ( PirlMul ( shade , ramp pw) ) ; 
CHECK{ PirlClamp (shade) ) ; " 
CHECK{ PirlConvolve3x3s (shade, 3c, )c+3)); 
CHECK ( 

VlPProjectOver (shade, reconow, 

((perm & WZ) ? slice cnt++' : {yoff2-slice cnt++) ) , 0, 
(WNORTH ((perm & WZ) 7 WY : WX) ) , 1.7 WRGBA) ) ; 
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pt__pw - scratchpw; scratchpw prevpw; prevpw = pt pw; 

pt_pw » grad_t; ^rad_t « grad_m; grad_m = descpw; destpw « pc_pw; 

CHBCK ( 

PirlCreateSubPW(&reconpw, &subrecon, xoffl, xoff2, yoffl, yoff2, 
OxF) ) ; 

CHECK{ PirlResi2e(e\ibrecon, .reconpw, 4, 4)); 
if (vol tbl->slice pos.z - vol_tbl->vol link->slice_pos . 2 < 0} 
if ({perm & WXY) " WXY) " ^ 

CHECK( PirlCircularShif t (reconpw, 0 , -50) ) ; } 
CHECK { PirlDeletePW {&subrecon) ) ; 
CHECK ( PirlResize*(reconpw, buf jpw, 4, 4) ) ; 

ROI WIDTH (vol tbl->roi> . ROI HEIGHT(vol tbl->roi) , number sic) ) ; 
if (perm & WXY) { , " 

xshift « ROI WIDTH { vol_tbl - >roi )- (vol_tbl->roi.y2 + vol tbl->roi .xl) ; 
xshift /= 2;^ xshift « (int) ( (double) xshift * sin(angTe)); 

FileHeader. width « PWXSI2E(buf pw) ; 
FileHeader. height = PWYSIZE (buT^pw) ; 
FileHeader .multichannel = 1; 

sprintf {filename, "*BV02d", f_name, fcnt++) ; 
fprintf {stderr, «5fs\n", filename) ; 
if ({fd « creat (filename, 0644)) < 0) 
{ printf {"%s cannot open file %s\n", 

FuncName, filename) ; 
} else { 

write (f d, &:PileHeader, sizeof (FileHeader) ); 
CHECK ( vxrr/copy (r>ut pw, tct, WRANGE, WRGBA, 
OxPCTJO, OxBFF, 
OxFCOO, OxBFF) ) ; 
close (fd) ; ^ 

} while { (angle += incr) , angle < range); 

CHECrK ( PirlEndLines { ) ) ; 

CHECK ( PirlDeletePW {&scratchpw) ) ; 

CHECK( PirlDeletePW(«egrad t) ) ; CHECK{ PirlDeletePW {&grad m) ) ^ 

CHBCK{ PirlDeletePW{&destpw) ) ; CHECK ( PirlDeletePW {&shade) ) ; 

CHECK ( PirlDeletePW { &reconpw) ) ; CHECK ( PirlDeletePW (tgradpw) ) ; • 

CHECK ( PirlDeletePW (&orevpw) ) ; CHECK ( PirlDeletePW (itemppw) ) ; 

CHECK( PirlDeletePW (&fauf _pw) ) ; 

return PIRL NO ERROR; 
error : return PirlEas'rerr ; 

coord3D_dbl *rotate_vec (light , light_dir, angle, perm) 

coord3D_dbl *light , *light_dir ; 

double angle; , . ,^ 

{ extern double cos() , asmO , sin() ; 

angle *- (-asind.) /SO.) ; . . v - 

fprintf (stderr, "light dir x %e v Ve z %e\n'', ^. 

light_dir->x, light_dir->y, . light dir->z 

switch (perm & WJCYZ) { 
case WYZ : 

light->y - light dir->y * cos(angle) - light dir->2 * sin(angle); 
light->z - light"dir->y * sin(angle) + light_dir->z * cos (angle); 
break ; 

case WXZ : . , . , , . ... 

light->x = light dir->x ♦ cos(angle) - light dir->2 ♦ sin(angle); 
light->2 « light"dir->x * sin{angle) + light:_dir->z * cos (angle) ; 
brea)c ; 

case WXY : 
default : 



noer ror : 
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light->x light_di3r->x * cos (angle) - lighc dir->Y * eiri(ahqle) - 

light->y = light_dir->x * sin (angle) + lightl^dir->y * cos (angle)'- 

fprintf (stderr, "light dir x Ve y %e z Ve\n« , light->x, lignt->y, liaht- 
retum light; ^ 
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It is understood that the exemplary methodology 
described herein and shown in the drawings represents only 
a presently preferred embodiment of the invention. As 
those skilled in the art will appreciate, the present 
5 invention is suitable for use in a variety of different 
applications, other than medical imaging. For example, the 
present invention may be utilized in radar imaging, machine 
recognition, and various other imaging applications. 

Indeed, various modifications and additions may be 
10 made to the described embodiment without departing from the 
spirit and scope of the invention. For example, various 
different shapes of structuring elements, other than those 
illustrated and described, may be utilized in either the 
morphological skeleton forming process or the 
15 reconstruction process. Additionally, various different 
criteria for defining the present membership of adjacent 
data points during reconstruction process are likewise 
suitable. Thus, these and other modifications and 
additions may be obvious to those skilled in the art and 
20 may be implemented to adapt to the present invention for 
use in a variety of different applications. 
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Analysis and Segmentation 
of 

EKgher Dimensional Data Sets 
with 
Fuzzy Operators 
for 

Representation and Visualization 

Wolfgang K Kraskt 
Nordurop Corporation 

Ahstract 

Hie segmentation and xepresentation of complex features in taigfaer dimensional data sets is 
of paramoimt impoctance for mnrhtM xecognition and hmman peiccpikgo^ image infanna- 
tion. A moltiresolntion and mnlti-^aspect data frprrwutarinn paradigm, £ moipbologicai 
skelttco, is used in tliis papex to provide a hierarducal framework for cffiriimr representa- 
tion and visoaliotion of data formacfaine recognition andbtman perception of datafeatnres. 
The n tiifratifm of fozzy Operators rTt^*^^^^*^* a basis within this tamewoxk for oganizing 
p^f^Tt Sy or fozzy sft s, of appioxlm ate infonnation by tnmiimiM» coverings or maximal sub- 
^ftKT- S-D 2r-D image data are used to demonstrate applications of these techniques on 
ht gh^ dimensional sets. 

f»rny«?ab> mathwnfl»traii fwnrphnifigy providcs an established basis and an algebra forfuzzy 
operators due to its representadonof fozzy nuq;soverasetsiipporLSpecificaUy»d^ 
tion of mocpbological operators in scale and orientation paradignis with tractable si^pport 
shapes provides an ordered basis for topological analysis and user perception of data. To 
eliminate predsioo loss grayscale morp^ogy utilizes CHily set operatitHis, requiring only 
con^MitBr addition and comparison. 

IntrodnctioD 

Thi% paraHigm ^ \^ T^n/wpholngi/^l «V^I»*ni « tntmriiir^H far hierarghigal orpnfraririn of 

fnzzy operators based on resolution. Operator smoothing and diOerencing opera^ 
aT*%#! ^ oitf^ frt ptipii n<V>^«1riw>lft f«it i ifff> r g pn»-cffntaringift mnilartp inner prodngTccpvoliv 

tional basis functions in the square mtegrable spaces. 

ry ot9 ective of this paper is to presen t an algorithm which preserves topok>gical xnorph^ 
and segments AMt^tt^ y ut t^igh resolution features to providp a 5ii'mp^**^f^^ T^py^ye^fo- 
rion far marfiiii^ wy^'r^itinn ami hrnnati viomitMrim 3, 21 . The fozzy sct algcbtaic tech- 
niques of grayscale mathematical morphology are tiie basis of this algoritfartL 
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Mathematical morphology offers a multiscak or multircsoUition approach for analyzing 
and segmenting higher dimensional daa f catmes aaoss a spectrum of neighborhood scales 
and lesoimions . Set theory is the underlying conc^ of mathematical morphology 

whirh p^i^m^ thg rtpreaentaHop of data with an open toiX)logy such that data anal ysis with 
an algf bfTi a n^ tppologrc^^ »^ « ™^ ™ contrast to direct analysis of data tx)ints . 
The extension to grayscale mathematical moiphology with firzzy set theory ^ faci ] i f af fi s 
the analysis of data represented as collections of convex functions over an open topological 
base. 

MathertRi introduced the algebraic set analysis of mathrmatiral morphology upon binary 
images to defme probabiHty distribotioDS of size and shape spectrarelat^ 
lie! sets eg stroctntingeleinents^. This approach to datarepr ^ 

sets rrtrTK^ff to rPTTT ^^^ft "^ »^ grgygcala or ftmsy mappings over sets. This approach devi- 
attt frrm thr rl antira^ of highiT rfmntictmal data acta perceived as aianrtom a nd / or 
discrete coUectioo of numbers or feature vectors to perform duster analysis upon 

For direa applications the algebra of mathnnariral morphology provides an expression of 
f i>amT f> riiTw#>«innality flndmorpholoyv for analysis and cl assifiratir o. In particular the tiac- 
lability of shape analysis allows the selective study of datafeamres of particular dimension- 
ality and size. Grayscale morphology extends this selectivity to l ocal cha rart r r istics with the 
quantification of texture and gradient as welLlhec^tability of soeasuring maximum and 
mxriimnm size extents is also possiblcby the diametric compariscM of shapes of stru^^ 
^iiTwifiHt at various sTflMft'^nrif"^*^* capability of measuring local intensity and a 
global morphology of imagery may be particularly beneficial for the quantitative analysis 
andmooitoring of physiok?gicalcfaangesintbepatfaologiesof turnorsinb ionirdtral ima gery 
andandtcpogra^ihicai changes for arid reconnaissance of evtroomental data. 

Thr i drfl f^thr sp^^rt^^ rfutnhntim of datarelative to abasis, such as theFocriertransform is 
as well suiyorted in the algebra of mat hrmariral morphology witii size and orienotion dis* 
tributicHis to describe images ^ ^ . T h ese morphological distribudoos re lat e shap e inf onnar 
tica to size aiKi orientation spectra of image data just as the Fourier traiisfom 
of sinusoidal functions to images " . Furthermore* the spectral distributions relate to U)e 
con^lexity of describing images with morphological shapes. Simple shapes such as mai>- 
mi ¥^- ^^^m^ parts typcally have shape distributions consisting of only a few spectra 
whereas natural biological structures ^pically have spectra over a significant range with 
frequency distributions of the spectra conforming to exponential functions. Hus function 
Tfr pfTHf f~^* f>b ^'^t^**^ Haf«A»ff HtfTv>tt«nn 10 of firactals ^ . The sh^ rharartfri stic of 
biolbgicalstnicturesispostulatBdasbotha&actalnature ^•'anda geometric nature in 
fact the Buffon Symposium of 1977^3 addressed both tite topics offiactal synthesis (Man- 
delbrot) and the morphological analysis (Sena) of biological stnictnres. In retrtiqxct it is 
tK^fnr^ p ifm«w <> fmwi tfaftitft mntnal intgrgsts In bi rima th rmatics thatmorphological analy- 
sis and fractal syntiiesis of biological sa u ct ur e s are interrelated by a common set theory. The 
cynryei* nf th\% d lsn w^ fitigutatift* the philosophy of a broad application and the unifying 
prospect of morphological analysis suggested m this paper. 

2.0 Methods 

The primary goal of mathematical morphology has omventionally been the representation 
of dag f catuTM by matrfaing if ti <>f r< *^ ^ i V «^ or imtarinn, the approach of this 
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paper is to r c p r rsCT t data features as deviadoos £ram sets. Tbe selection of the apprpphate 
diznessioDality and shape of a structming element or kernel set for mcHphological analysis is 
fundamental to ttc techniques ctf this jtaper. A rigorous develofment in henceforth devel* 
oped from topological coostnicts to operator analysis and finally dxmetttional filtering. 

2.1 Topology and Structuring Elements 

Tbe concepts ctf topology are fundamental to tbe analysis of image data relative to varieties 
of sets disdngnished by *>wp^ orientation* size and position. Tbe topology used is founded 
by the inm or fimtiliesofsemiixxms on vector q»oes. Tbe text ofFoll^^ 
pensable rr f f rr n r g in tbe development of this formalism . 

2.U Definition : A vector space X, over a field K, is closed for 
tbe following coniinnous linear m^ipings given any 
X€K, x,y €X: 

i) vector addition X x X -* X, (x,y) — x + y. Vx,y€X. 
H) scalar mnltq>licati£Hi KxX-^X« (X,y)-*Xy. 

Tft Hf^gin, m p ^f^ ior, w will Ki> ivwi tn anal yre images based on neiyhborfaood pro p erti es 
over a vector Space X. Le. 

2.1^ Definftftoa : A topology on X is a family ST of subsets of X 
with the following prcpeities : 

i) X 0 € r 

ii) if the cardinality of A, tAI S K « INI , 

and (U.).c^ C r dien U Ua e g^. 

•ex 

iii) if lAJ < K and {U.}.«a C 3* then H U. € S', 

nou bene : lAJ denotes the cardinality of a set A, 

N denotes the natural numbers, |N| - K. 

Hence 9* on X defines an topological sfiax (7*. X). This paper will consider topological 
spaces with countable bases, defining a separable space. 

2.13 DeflnltioD : A neigjiborhood base. JTCE) of £ C X 
has tbe following p rop er ties : 
given a set T(E) « {A : E C A. A e T), 

i) JTCE) C rCE) c 

a) V U € TXE) 3 V e Jr(E) such that V C U. 

Definliioa : A base of a topological space OC^ is a family 
of sets 8 C y such that V E C X, fCE) C 8 and for any 
ndgfabcvbood base JT. of E, X(p> C 8. 

Tbe vector q»oe owr this algebra wiU be tbe closure of X. This pqxr win only a 
closed vector spaces, denoted as X. 



wo 97/09690 ^ " PCT/US96/14500 

3^ 

AnalyiH and ScgnntBtloa with Fuzzy Operator* for Repre»«irt*tk«/ViK»lix«lion 

2.1^ Definition : Given a topological space. (X,90, tte dosuie 
of ?* is a o - algcbrm. JL, such Utat : 

ii) lAI s K. and {U.).eA C ^ then (J U, S ^. 

iii) lAI S X. and {U.}.eA C ^ then H e a. 

(note : If X • X then X is complete relative to r, a Banach space) 

To simplify analysis of itaages, restrictions on countable bases a of 5* mnstte conadenri. 
Fortbeimagedaiansedin this paperXisanEodideanveaors^ either R'or T^^R^/ 
aZ».ftortaiesi«o. iftopotogyisconth«oas,orZjorZ^/nZ3if topotogytodisar^ 
iMiesi»e7€ z. Tte discrete topologies onZ'orZS/nZ^haveaconnailebasecrfpomts 

sm^ets with finite collections of potots which defines a separable space. For a connnnons 
space it is necessaiy to know if a base can be developed from a countable set of scales and 

txanslatioasofanopenset,E.widianit xadios, RCE) = 1. see below. 

XW DefinitioD : A Msoinorni on a vector ^nce is a 

mapping X [0, «) with the foUowing properties : 

i) p(x + y) s P(x) + P(y)- V X, y e X. 

ii) pOx) s PUp(x), VXeK,Vx€X. 

if Oii) p(x) - 0 iff X - 0 theo p is a norm . 

denote d : X 10, •) as the EocUdcan norm. 

2.1.7 Definition : A diameter* of a set E C X is 
d(E) « $op{d(x - y) : x,y e E), 

2.1^ Definition : A radioa, of a set E C X is 

ROO - inf{x e E : sap{d{x - y) : y € E}) . 

Anther set property is now convenient to introdace, the concepts of open sets. Topological 
spaces in particnlar are construOions of open sets. 

2.1-9 Definition : An open bafl is a set B^x) such that 
B/x) - {y : x.y € X d(x.y) < r. r € R. r > 0). 
vA^at d is the Eudidean metric B,(x) C X. 

ZIM D^inMon : The interior of a set. A, is the union of opm balls 
contained in A, 
A - U(*^X:Bg«(x)CA, for some 6(x) > 0). 

0 

2.iai DellnWon : An open •et, U. is a set such that U - U. 
Given the definition of an open set ainore gencxal concept is the definitioo of aneighbor- 
bood. 
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2.1.12 Definition : A Deighborhood, A, of E C X, is a set 

o 

siidi thai £ C A. 

With open sets it is now possible to rtrfrnr mmpRa sets, wiiich will ccmplcte the intioduc- 
tiCQ of primary topological concepts. 

2.1.13 Definition (Heine — Borel) : K C X is compact implies thai 
given K C U ^» • where {Ual.eA ^ an open coUectiozi, 
3 D C A, IDi < K, soch that K C U U. . 

2.1.14 Propoaition : If K C X is con^Mct then K is bounded . 

Le. K C BJixX for scHne r ^ 0, x e X. 

2.1.15 Propoaitkio : If E C R' is bounded and dosed tlien 
E is cnmpart, 

In topological analysis it is convenient to have opca sets with an additional bounding proper- 
ty shared by minpa r t sets in vector s p ac e s. 

2.1^6 Deflnidcm : If E is compaa and E C X is open thai E is 
precompact hence, Br«(x) C E C B^x), 
for same r.r' > 0, x € X. such tliai r' ^ R(E) s r. 

Since tiie analysis of this paper is restricted to countable bases of X, be it or , the topo- 
logical properti es generated by the translates and scales of a preconqsact set £ are estab- 
lished. 

2.1 JL7 Pn^MMitioo : A simple base of a topological space 
CX,50* is a base 8 C if such that 3 E C X, 
E preccanpact, 8b - (£.(§) : a e Q, p e A}, 
E,0) « o£ + Q is a countable subset of R, 
A is a countable aibset of X. 

Proof, assume 0 € E otherwise choose x € £ and 
consider E' » E - x which contains (0) , also assume 

R(E) - 1 otherwise set E' - rJqE. 

Given UC?", VxeU3 B,(x) C U and for scnne 
r > 0. E^>(Wx)) € 8 such that "E^yi^x)) C B/x), 

Le. d(x - P(x)) ^ J and a(x) ^ J, 

then U - U E^(Wx)). hence 8e is a base of 7. 

sCU 

The condidon for a generating set E of a sim|de base 8e to be pircnnipart can further be 
relaxed to i nrl w ^ any set assuming the second condition of the Baixe category tiieoxem is 
cari^^^ h#mr^ this tbeorcm is stated below. Fust the defmidon of a dense set is required. 
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XiaS Dcfinitioa : If A C X . then A b dense tdativc to X 
if 7C - X . 

2.1*19 Baire Category Theorem : Let X be a complete metric qace. 
i) If (U.l.eN is a countable coUecticm of open 
dense subsets in X then Q U. is dense in X. 

U) if X « U E« then for some k G N. is a neighborlKX)d 

of some point in X. is the closure of . 

These generating sets of the simple base of atopology arc commonly xcfcr^ 
ing elements inmathemaiicalmQrphology. To reduce complexity of amlysis these structur- 
ing eloncnts are further lestiicted to be convex as well as prccon?>acL 

Xl^ Definition : E C X is convex then V x, y S E. 

Xx + (l-X)y€E» VXe[0.11. 
GivenapreccmqHcta)nvcxsetE$cvcralprc)pcxticsusefidforfu^ 
tional analysis will be characterized. The following Tn^irhrmatiral formalism for convex 
function spaces has been extended from that established by Conway 3 . 

2.1^1 Proposition : given a prccompact - convex set E 
3 a tmique point c € E from wtdch tiie radius of 
E can be detinofi, E C Bb(s)(<1) * 
Destoce this point c as the center of Radius of E. 

Pnxrf. consider R,(E) » $up{d(x - y) : y € El. 
since R(© - inf{R.(E) : x € E) a x e E 
such that R,CE) » R(E). 

Suppose R^ « R,(B) - R(E) with q l . q.t € E. 
Hence V x € E, dt(x). d.(x) < R(E), where d<x) - d(t - x). 
Hiis implies that £ C BK(q)nBR(t) , ^ ^ 
But E C B.(q)riB.(t) C B^OO . X = ^^ . 

X € E since E is convex and R(E) ^ R,(E) < R(E) 
with q t, a contradiction. Therefore q "> t c 

Furthermore ifEishasasufBdcntnumberof axes of symm«ry die unique center o^ 
a member of the precompact strocniring clement E, 

Proposition : if a prccompact - convex set E has hyperplanes 
of symmetry (Piir.i with a unique point of inicrsectica. c 
this pcrint is both in E and center of radius of E, R,(E) « R(E). 
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Proof, given P € {Piir.i , By symmetry property if 
X G E then 3 x' € £ sudb that x' is the reflection 

of X across P. ix. X « l2-L£> e P. 

T = Q^x) = inf{d(x - y) : y e P} . 

Hence if c is center of radios of E, R(E) ^ KSE\ 

then the reflectioo of c setoss P, c' is also a center 

of radius of E R..(E) = Rc(E). by the previoos 

propositica 2,1^1 c is muque and c » c\ by symmetry 

this c 6 P . But P € (Pjr.i is aibiaary 

hence c» O fL,Pi a nniqne poinL 

Now choose x € E with x ^ c this ixiq)lies that 

X ^ P e {PjIH, , firam above x', and Qr(x) e E. 

By induction c = (O**-* ^i)^*^ ^ ^' 
The stmctnring element is selected to resemble the mmpart neighborhood form of the data 
struct u re s of intexesL Rv fractal data ttiis is particularly convenient since the same primitive 
form and structure may be repeated over a broad range of scales. Hence the sdected strucnir- 
ing elanent f onus will locally resemble the features of interest in compact neighbcKhoods . 
Therefore, a coDpaa 3-4> ^jhere and 1-D arbitrarily oriented line segm 
Ki /wM>rf;/ni tmagi*r y jttrngtnring elgments to represent the neiithbochood of brain smtace 
structures vascular curvilinear stractnres respectively. The data analysis of com^dex 
tmag #> f j^fltttr ^c with oorrrpart anri rrm year stnictiiring elements is casilv implemented on t>ar- 
allel processing digital ocRupoters. 

To demonstrate the shsqse of the H t ]«f *i t'*'imat/> sphrriral structuring element* note that the 
structuring m Figure 1 is eight sided and symmetrical in two sets of four directions. 

Note thaf the of symmetry of this structuring element require the center of radius to be 
within the structuring elementp propositicD 2.122. This sh2^ possesses a quality of direc- 
tional isocropy which is ideal for size aiui on vexity analysis of 3-4> image feattires h^ 
itislikewisenotdisposedtoprovideananalysisoftheorientatioo ofxmage features. Further 
analysis of the curvilinear features of 1-Dfonns, such as the stexiesaiid veins of 3-4!>MRA 
ft/v jiiieirinwe re quiics Gfiented structuring elements resembling 1-D line segments such as 
elongated '»"ip«^ Syumietry about a central point is common property of all structmiug ele- 
ments used in this paper and thus the rrK^Wv^*^ system of the structuring element originates 
£ram this cental point, or ceiuer of radius. 

Stnrr imagf^OTfti ^"**'**Tr ^"'^^ '*^^^ 

rt/wl rywjiwt tf*f> tM*w> ff wmaKqn< nf ttwatynic mntf extended tO the rTpfrVTltatinn Of 

functions miqsped over X. In discrete topologies with a countable base however, a sphere 
ftfvt mhs^ r^tinn fw yt famnru niat in a digital repr rOTitariff¥ i Het>ce thc Structuring 

f>ii»mf«f ap piff^ in this paper digitally ap pr oxim ates a sphere 

A convex fnncticn» will be introduced in the next sect i on , m^iped over the structuring 
^u»Tfyfi» domain lo extend the topological basis to a functional basis. If g is noiHCOQStant 

tfaf^ atnagniring ^l^tmn f \ yfflin4yri<^^ are firrrififri , prnnttlrn g more 

1**^^ f n * ^' ■ ■■ w%i>tin>i r>f ft «ph^ flnrf pr^t^ n^iihi\ity With yencralixation m d a t a a nal vsis. 
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Several generalizations of the set basis definitions are now presented for positive function 
m^jpings over the basis set, primarily, the stracturing element diameter and radius will be 
that of the support. 

2.1^ definition : the diameter of a positive mapping 
g is d(g) = sup {d(x - y) : g(x), g(y) > 0 }. 
Assume d(g) - 2 for g tmless otherwise stated. 

2AJA dcfinitiDn : ttie radios of a positive mapping g is, 
R(g) - inf{sop {d(x - y) : g(y) > 0 ) : g(x) > 0}. 
Assome R(g) « 1 f or g unless otherwise stated. 

2^ Foadamental Algebraic Operatiom 

Tliere are two ymiary neighborhood operaticHis of maf hrmarical morphology, which are re- 
latedby set complementation and are critical for the analysis and development of tiie mor- 
phological representation. These c^wations are dilation and erosion. Conveniently theho > 
mogeneoos cascade of these neighborhood operations increases the scale of the structuring 

element without precisiim loss. To perfonn these critical operadons of dil^^ 

the stinctnring ciememinast be translated over the nnage domain. To fomulai^ 

tionsadcfimtion of the ftmctitKis and the bases winch siqjportd^ 

must be derived. Tlie approadb of this developmem is to extend the dcvelo^am^ 

section to function bases which provide analysis of function or image mq?pings m higher 

dimensional vector spaces X. Dilation is defined witii properti 

operation of erosion is intro duc ed 

ZXl Definition : g is a convex mapping over E, £ € Be , 
wiOB E is a aaiq»ct convex generating set of B^. 

a, - {gx + P : P € Q. A € Be) 
For ease of notation ag » g^ , 
Le,ag(ix) -g^x) V x e E. 

hcnce» ag(aO) - g^O) « %M V e > 0. 

With the basis, Sg, tlie objective is to approximate functions to arbitrarily precision in the 
space of ftmctkais witfi Ixbesque measure over the topological vector space 

2JL2 If/:X — K where K-R,Z, 

Consider CV- {x:xeX l/J > a), fimhermore 

since/is a function if a > b, C/C C/, 

then Irt. - inlla 0:m{C-j) - O), 

alternately. ^ ^ 0:m{C*;) > O), 

m is die Lebesqoe measure on (X,^). if/€ L* then I^. < 
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Docc : to sin^slify notadoD assume (/] = , 

A more gemal class of functions can be admitted if these functions are considered inte- 
grate over precosnpact ndghbcalioods. 

? definition : / G L£c implies V pre comp act U C X« 

-= infja ^ 0 : m(Ci,;) = 0. U C U' e ST) , 

wbcxe q^/* u'ncv- 

Withoot loss of generality/can be xestzicted to a positive funntiop over X. It can now be 
< hftom that tfaf ^ bay. P^E* <^ ^rrmrimflte a function, /, arfaitrarily closely in the nonn, 
Theiesolts win then be shown to extend to tiszng the function base, Sg. 

20-4 given F C X, the measure of F is : 

m(F) « inf{m(U) : F C U e iT} . 

Using the L~ nonn there exists a mapping of L*"ioc into the set of upper semi-continuous 
functions, ^ , with function bases, « supported by the tiie precompact topological bases. 
8e, of section 2.1. 

20-5 definition : If /€ ^ /:X-^K, then 

•.a)) is open in X V a € K. 

To faHKtatff further development a new functicxi will be drfinrd 

20.6 definitloo : A " + «*) X^g • Xb is the 

characteristic function : 

f 1 . if X SB ^ ^ 
X«W ^|o.ifx«£B' forBCX. 

The most fundamental conation for the image analysis developed in this paper is now 
introduced, the dilation nuking relative to a structuring element E. 

20.7 if^n^irt^ : A dilatioii mapping, T, perfonns 
the following operadoo : 

Tim - tra - ^ to - to - ^- <i) 

E.(x) is the reflection of E,(x) about x. 

likewise ^ is the reflected mapping of g, about x 

with E«(x) the sappcxt of g^. 

note : the analysis of this p^per will use symmetric g, 
soch that, s' - g: . 

Denote T\ as T. if the structuring clement g is understood. 

Th ft introdiictiaD of dte function su^sping of dilation permits the proof of several interesting 
results for image sialysis : extcodibility, continuity skI onler preservation of 
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2,7* proposition : 

i) TJ r^if E a e' . 

ii) TJ (x) is contiDUoas with rKpea we, 

V e > 0, X te X. 

proof : 

(i) ^(y) a fc.(y) V y e X. if e a 6' . 
since g is omvex, 

r^(x) « Ira - e . where g^(0) - e , 

A.(x) - E.(x) \ E..(x). E e'. 

Therefore. T^a T../ since FJ'Cx) r../W V x S X. 



Oi) r*/ is coDtiimous from bellow : 
snp{r»/(x) + b : b < e) 

»= sop sup [a a 0 : xn(C*y^) > o) . b < ej 

- sopfa fc 0:ni{C*/^ > O} . since^ for e > b . 

- r,/(x) + E . 

" ' r./(x) - snp|rwr(x) : b < e) . 

r,/ is continuoQs firom above : 
infjrfcrW + b : b > e) 

- infjinffa & 0 : m{F.(Cj;;)) = o] . b > eJ 

- nif|a 0:m(c|5sj/i) - • 

Le. - n . E»(x) - supp(/i) . 

- inf|a 0:m(ci53/i) - 0 j 

- l/^(x)iB3jj - 8iip(lr(x)|iw.o . 

- ^(x)!^ « r. /(X) + E . 

or , 

r,/(x) - inf(r^(x): b > e) . 

ny rf^wirinn tho. ha.<M ftmetioM of a, are boonded coovex mappings ovg a preco mp a a 
convex support With dcfiniiions of a nonn it is now posable to develop a nana with req^ 
to the convex ndghbortiood properties of the snppcct of g. Le. E s siipp g is convex. 
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dennltion : A set £ is halarrfl if Px e E, whenever x € E, 
and IPI £ 1. for p € K . 

A >«i^"^^ set t^ T^*^* symmetry of tt>e set aboot any byperplane through the origin. Hencs a 
precompact-coDvex stracmring clement set which is halamnrd comains its center of radius, 
which is as well a center of symmetry, as specified by proposidon 2.122, 

22J10 imposition : Given a preccmpact, convex and 
balanced set E C X , {0} G E . then 

Pb(x) " inf{X : x € X£} , is a nonn on X. 
pnxif. Gtven a precompact convex set, E C X the three 
properties of die nonn are to be demonstrated. 

i) Pe{x) -= inf{X >0:xeWE}2sOVxeX since 
3 5 > a 0 e B*(0) C E. and ^jjjpf-^x € B^iOX 

or X e ^^^> ^ for any £ > 0 and V X e X. 
o 

tbexefcre 0 S p.(x) S V x e X and Pb(x) = 0 iff x - 0. 

ii) Pb(Px) - mf{X > 0 : Px S XE) 

-ipimf[^>0:xe^E|-iPli%(x) 

V p e K - {0). 

Ifp-O, px-Oand Pe(Px) -= 0 by 0). 

nJi. that if p 2 0 . tiie set condition can be relaxed. 

iu) Given x,y SX e > 0. 

''-S5^-'^=S5^ sEsi«.prfxo.,.^')<i. 

Using tbc convexity of E widi t) E [0, 1], 

then tpc' + (1 - Ti)y' € E, and p^dix' + (1 - rOv') < 1 . 

Now suppose 

n « pg(x) 4- c . 

^ Pe(x) + PE(y) + 2e "™ ^ ^ 

I>e(x + y) < PbW + PB(y) + 2«» bot e is arbitrary, 
bencc p^x + y) < p«(x) + pe(y) V x,y € X- 

Thisnoan in X, proposition 2^10, can iiow be extended to coiivcx mappings 
compact-convex-balanced sets. 



7 7 |t Thmran : Given a decreasing convex mapping, h, over Pe(x), 
h: [0.11 — 10.1] . h(0) - 1 . the composition g -hQpeisa 
convex mapping over H, with die pioper^ . 

^ ' ^ X 6 £E. y € tiE. 
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proof. Given a precompart - convex set E C X. 

" zXx\ , 11' •= 1 - e'. x' «=• y' " e + n • 

^ e'h(£^) + ^ bCPnU') + PB(y')). 

sinee b is convex. 

By pmp. 2.2.11 ^x') + peCyO & Ph(x' + /). 
f^KitiiKi with the decxeasing property of h . 

h(p.(x') + piKy')) « hWx' + yO) = ^{tt^ ■ 

hence, egg) + tig(^) s (e + »l)«(fT^) • 

This prapeny of g iiiq>lies g is convex. 
Since Pb:E— [0,11 and h : [0, 1] — [0, 1] . 
titen g - hOpB:E— 10.13 • 

2,7 coroIIat7 : Given g a convex mapping over E 
r,. g. + e' - r, g. + e - g.*,. . 

ISDof. Since g is coavex g is ctmtinuous, hence 

r.. g.(x) + e' - si?>{e' g{^) + e : y « E.-) 

- (e' + e)g(7-J^) - g.*. (x). by theorem. 
TheRfoTB r,. g,(x) + e' - g,*.-W • 

2J2.13 : Given e' — e - 6 > 0 and / G LHU , 

then r./(y) S r./(y) + gl(y) s r,./(x) + 6. V y e E,(x). 

proof. If /€ Lta and e' > e then r,/(x) S r../(x) , 

by (i) prop. 2JL8 . 
Sete'-e>«&>OaDd choose any y € Et(x) 
«b« + fcCy) s g'. by i»op. 2ill. hence + &(y) a; yj, . 
wUch impUes |^ -4- g^(y) £ ICj and 

r./(y) + gI(y)-6 s r,./(x) . 

T1ico««m : F, : L£. — «U C L£, V e > 0. 

pnoL If r./e«Ufor/eLi then 
AT./- {x € X : r,/(x) < a) is apea m X. 
Owoae X e A* rj' then r./(x) < a: 
Choose a' > e then T,f(y) s T,.f(x) + 6 . 
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for any y e E^Cx) and e' - c = 6 > 0 , by (i) lemma 2.2,13 . 

Choose 6 small enough such thai r../(x) + 5 < a, 

and hence FJiy) < a. V y € E^{x) C A» . 

Since x £ A' F./ was chosen aititxarily this 

implies that A* T,/ is open and F,/ G 

Furthennore since / G L£e is aibitiaiy it is shown 

F, : L£, — HI C L£e V £ > 0. 

The now fiannally introduced image transf onnation of dilation not only projects the function 
space L«« into the upper semicontznooos functions but xestricts the neighborhood m^^piag 
properties of tliese functions. 

20.15 Definition: If/ G at then limF»/«/. 

2^16 Lemma : Given/e «U . c > 0 . then V x € X , 

F./(x) « sup{/S(y) : y e E.(x)} - t . 

piDof. Suppose /e *U then by prop, 2JL13 V x e X , 

F^(x) + 6 sup{ F..4/(y) : y ^ E»(x)} . 

hence F^{x) + e $op{ /(y) : y e E.(x)} . 
But since /J(y) - c ^ /(y) . 

F./(x) ^ sup{/S(y) - £ : y e E.(x)} 

St i;a - « - rj^w . 

hence F./(x) - $up{/;(y) - £ : y e E.(x)} . 
2,2 17 Lemma : Givcn/€ Li . £ > 0 . then 
F./- F./, where/- hmF./€«U . 

proof, (a) If / e Lte tbeo by definition 
F./fe/ + g^-rc on E,(x) a.e. 

where ax. means almost everywisere except on a subset ol E«(x) 
oi measure (definidon 2JIA) zoo . Hence 
/(x) « limF./ ^ lim (y) + g'(y) - £. y e E.(x)} ax. , 

therefore/^ /ax. and F,/ 2 F./ for aU a > 0. 

(b) Now suppose x € X such that F./(x) > F./ (x) , 
by lemma 2^16 ,3 y € E.(x) such thai 

F./(x) fc/^(y)-£ «/Cy) + gXy)-5 > F./(x). 

but by kmma 2il3 /Cy) + gKy) ^ F./ (x) + e , 
a contradictiODv 

Iberefore (a) combined with (b) tn^lies F./ (x) - F«/ (x) . 
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2^18 proposition : given the dilation m^ings T,, T,. , 
constructed with the convex structuring element g , 
the composition map has the property : 

r..r./ -r,.*./ v/€Li. (2) 

proof. Suppose /e LS, then by lemma 2,2,17, V x e X . 
r,/(y) - e - $ii>{ 7J(0 : t e E,(0) . y e E,.(x). 
T,TJix) + e' + e 

- sup(sup{ yjft) : t e E,(y)) + fc.(y) : y e E..(x)) 

- snp j fit) + sup(^(t) + g*.(y) : y e E,.(x) } : t € E.(y) } 

- sup(r(0 + gi'^fi) : t € E.(y) ) 
-r,.*./(x) + (e' + e). 

HwrefcHe r,.r,/(x) - T,.*, /(x) V x € X . 

Wth the now standing dUaxion <^>eiaior properties on L~ioe 

continooos functians of *U a new set of fimctions can he derived. This dass functions is 
called the t^pttrirted lower semi-^ontinooos fiactions of X-Tliegen OT of this d ual fu nc- 
ti(a dass is now introdnced by fl»e function transfonn operation of erosion, defined &om the 

opoation of ^latim. 

2,2J19 Dcflnltica : If /feO, /GLtthen -/SLSe, and the 
opetation of erosion is : f , /(x) - - T, (- /Xx). 

The operation of oosion is dual to tliat of diladm and introduces many additional possibili- 
ties for morphological data analysis. 

PropositioD :If/sO, /€Ltthenf,/€X.Ve>0. 

proof. Suppose/SLi, then -/€Liand 

r, (-/) € «U V 8 > 0. 

Tbeiefoie, AT, (-/) is open and CT. (-/) 

is dosed . Reverang polarity again then shows 

(x € X : - r, (-/Xx) > - a) is cpen and € A . 

S ince dilatiaa and erosion aie duals, each property of dilation also has a dnal property with 
erosion, which can be proved with this duality principle, hence proofs win be omitted. 

DcfinhkHi: If/€ I. /2 0, lmf./-/€i. 

'>'>'y^ imposition : 

I) tjS r,/if e e' . 

ii) tj (x) is continuous with xespea tot,Ve>0, x€X. 
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7223 proposition : given the erosion mailings F,. F,., constructed 
with the convex structuring element g, the composition 
has the property : 

f..r./ -f..../ V/GLi. (4) 

If the finite values of X and g are binary valued these definitions of erosion and dilation cor- 
respond to dseir binary comterparts. 

With the operations of dilation and exosicm the less destructive operadons of opening and 
closing may be performed. These operations of dosing and openmg conform, respectively, 
to tfa g definidons of ttiOTmtTm endostse and Tn*«iimiin cover of a dataxn^pping for unions of 
the snuconin g element as it is translated. Opening is of particular xnqxzrtance in the develop- 
ment of the moipliological ^v^^^<^, gxninxmal image rcjnesentation, which is ^fi^rd f r om 
deviations of the image frimi the sh^)e of the structurxxig element aat>ss a spectrimi of so 

2224 Definition : fiizzy Opening, (5) 
« F.f,/, g is pinuj npa n — convex, e > 0. 

properties : anti-extensive C / . 
idempotent 

increasing / C * C . 

2^2S Definition : fiizzy closing, (6) 
f, F./ , g is precompact - convex, e > 0. 

properties : extensive • 

idempotent ^ * 
increasing /C/i-Z-c*-. 

The strucniring elements of this paper are limited to arbitrary transpose invariant convex 
shapes which are selected to measnre the size or orientations of image features as erosion 
aixl dilatiODS are successively cascaded. An interesting property of dilatiCHi and erosion is 
that a bomogeneoQS successiont>f these respective operations is equivalent to one operation 
with the oanposite dilation of the entire gn>up of structuring dements, eq. (2), (4).!^ the 

m wiiwwtrH ftf o rr»tw>y r^g^ U ti^^ft^ Ky ft gnr/y«tnn nf gfnftinng with a spheiical 

struonring fV****^ until the object representation is deleted. Likewise the maTimrrrri diame- 
ter of a convex obj eamay be ascertained with a soccession of dilations with the same spheri- 
cal strncttning elemenc 

Impleinentationofmorphological opening operatioiis for diazxietric analysis h^ pro- 
vided asegmeiitaticsn tool for the unobscured visualization of in tracraxi^ 
nal organs firom X-ray CT and magn^tir^ resonance imagery . In paniailar neurolog- 

ical surface anatomy is seginented from 3-1) acqtitsitions due topnsperties^ con- 
vexity acn»s a spectrum of scales relative to the concave 2-^ surface feanires of the 
oping extracranial tissues (Hgnre 3). For grayscale morphology, convexity also extends to 

t*w% hfi ght>>t 5farirrft withm f <rnif pf pft^fffjtf^ Siooe cxosion Bud dilation Or Opening 
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and closing are lelated by set comi^enientatiOQ it is equally possible to segment dark fea- 
tures by use of dilatico and closing (^)eratioas. This property is exploited to segment the ven- 
tricular systooa of the brain (Figure 9). For MRAsnidies the surface stn of 
S-I> erosim is crploitcd to ^iw«»natff extracranial c^illaries for unobscured visuaiiza- 
tioD of the intracranial arterial and venous tree. Unfortunately, most plications of mathe- 
matical morphology, i"iiriv#> the methods introduced xq this paper, do not extend beyond this 
level of morphological processing and hence destroy most of the high resolution detail of the 
segmented objects. 

Z3 Morphological Skeleton 

Thf Tqfff ^ti it^f" "f fiflp-cmvgx or concave tryological characteristics, such as curvili- 
near features, over a spectrum of neighborixxxl scales is supported in mathematical 
morphology by a fV^ii'tn^ryflriftn technique. The morphological slrrlrtnn representation of 
an «mag i> tie rVwJnpftd hy the coUcction of successive shape deviations of an image set from a 
stmrnmng glgmcnt aT^*^ ^ «pi>rTnn" cr-aV-c Pnr arteries and veins this maps a uniform 
diametric representadoo iR^ien the circular structuring element is used. 

It is apparent £ram the exanqples that morphological c^Kning and closing conations destroy 
features, withsnanimage, which deviatefromti&e shape, stzecx-orientationofthe structuring 
element To preserve these complementary framres the morpholopcal skeleton is used. 
Thus image features which deviate from the properties of size, shape and orientation of the 
structuring element are preserved witliin die morphological skeleton. The construcdon of 
the morphological skeleton of an image proceeds by <fi£taencing successively laiger scales 
of morphological erosions until the image erodes to a monotone intensity, or successively 
larger scales of morphological dilations of an image firom a largest scale, which just dilates 
the image to a monotone intensity. 

Thus y^^rm^A f eauffes widiis a neighborixxxl scale (white) as well as minimal feamres 
within a neighborixxxl scale (black) are rep r ese nted within the motphological skeletcn. 

To represent the white ^t^i^tnn ihe anti-extensive operations of erosion and opening are per- 
formed. 



23A DellBhion : SKfr f> » (T./) \ (f./)^ . 



(7) 



S*(gJ) - 21 ' skeleton. 

Alternately for a b^ck tV^i^^^n the extensive c^serations of dilation and closing are per- 
fonned. 



23a Definition: SL.(g,/) - (r,/ )''\ ( T./ ) 

lim inf 

b- {o: infers - supV^l. 

F<r digital faster foanats of i»xels or voxels e is eqnal to 1. 



(8) 
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S-(g, /) = SLuCg. /) . is the black skeleton. 

Z33 Definition : S'(g, /) « ^ S^Cg. /) , (9) 

the total skeleton. 

An important feature of the skdeton is that the original image niay be reconstructed without 
precision loss with dilations and additions of the positive moipb(dogical skeleton compo- 
nents SaOO or with erosions and additions of of the negative skeleton components S^oOQ. 
For grayscale images tbr^ equations are expressed in arecursive form : 

23^ Definition : « TJJ^^.k + S;(g./) . / -/o-, . (^O) 

white skriercmal reconstruction. 

Likewise to represent closings of the image the negative components of tiie skeleton are 
used. 

23JS Definition : = t^f^^'^ " S!.^(g. J) . f » . O^) 

black skrlrrnnal reconstructioiL 



2.4 Fnzzy Connectivity 

Connectivity isacharacteristic of tbemorphok>gical skeleton preserved relative to the corre- 
spoiKling image structures if a convex structuring element with finite surface ctzrvatures is 
sfilmcd. The ca^sacity of the morphological skeleton to preserve ti« connectivity of image 
stractures is exploited by the algorithms of this paper. Fes- digital imagery and fuzzy sets the 
a^niai definition of continuity, however, mnst be defined. The continuity of a skeleton rcpt^ 
scntation is definedrelative to the str u ct u r in g elements used to make the morphological skel- 
eton. To faHiitttfi* this definition the introduction of a modified Hausdorff metric provides 
the Tngan^ of measuring tiie connectivity of image structures. The Hausdorff metric was de- 
vised to iiieastire the distance axid dissixxiilarity of sets, anoiiliiiear set aiiak^ 
cmelation. H&is metric is precisiely drfineri by mathematical mcvphology if a circular struc- 
turing element is utilized. A modified metric may likewise be devised with variously shaped 
structuring elements of relative convex and syxxunetric ctiara ct e r istics. Proposition 2.2.10. 
Thu metric likewise generalizes to a metric tft*'^"*^ between functions with proposition 
2ill. 

modified Hausdorff metric between fimctional maps/and h in restricted ^ • • 

2^1 DeAnitkm : Modified Hausdorff function metric (12) 
AJif.h)^ inf{a : sup(/A A) :s inf( r./A r.A) ) . 

if/is connected by distance e > 0, and elernent g then dg(/(x).yiy)) ^ e for .aL: points 

Fgt tiie 3-D digital image applications of tiiis paper ccsnnectivity is limitrrt to a tmity radial 
v^y^ of tite 3x3x3 ba«* structuring elements relative to tlie Hausdorff metric. 
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2.4«2 Proposition : 

if an image/is distance e and element g connected then : 

So( C/. (13) 

proc^. So(g. /) > /u.* are distance e and element g connected , 
by definition 2.3.1. Likewise S,(g, /), ( f ./)^ (14) 
are distance e and element g connected. 

This pnqx>sition is e^qiloited in morphological reccmstnxcticHi techniques based on fuzzy 
connectivity. 

Thr ri'^ir^^^CKH^ ein^utrm « fli cn g rgriiiced T rxTTPVintatio n of tiie image, since due to prtrp^ 
enies of opening no set intetiots exist, which is ideal for image compcession ^ ' , and analysis 
of multiscale curvilinear featnres of images, indep en dently from diametric f eatnres. For 3— D 
fig^ a sets ti&e morphological y^f^^^^^n consists of fonns which are immeasorable to a 3-D 
stnictming element snch as 2-D ^iirfar?*. 1-D curves and 0-D point noise. For curvilinear 
fiMitTiw <i* £trimtariq n t bT^ftr phningical skeleton IS au ideal representation since arteries and 
veins are rep re s ented as 1+D corvilinear fonns widi nnifdnn diameters of thifimrss immea- 
sprablc to the scales of thg stroctnring r^mmt nwl Th« is dug tn the rthnmsional filtering 
rhar^/^i^TTft*^ of Opening f eamres of lesser r<tTnii>gtnnai measure are consistently cov- 
ered by collections of structuring elements of larger eudidian dimensicaality and thus fil- 
tered (Hgure 5) andplacedinto the morphological skeleton. The ability of tiie morphological 
yiffii*^f wi to rep r es ent features over arange of neighbodiood scales and within a range of inte- 
ger ritmmctftnai bouuds permits the representation of fractal features or fractional dimen- 
sion, wiiich is demonstrated in the results section. 

Z5 S^mentntioD of S-D Features Based on Fuzzy ConnectiTity 

Tnt rarTaniai fcamrftit of tf^ tgam flffK gegmgntcd fmm MRI argni^tioD based upon a morpbo- 
logical filter of their large extent followed by a morphological reccostruction based 
upon fuzzy coimectivity of smaller scale features from the grayscale skeleton components. 
This provides the restoration of the high resolution features, particularly brain surface and 
arterial features. The size filter is based upon ravartrd erosions, eq. (4), with a spherical 
gtrn/fwrmg ^itmmt Figtm tmril thg ertragranial tissnes are deleted due to the lack of con- 
vcjutyand diametric properties of the scaled structuring eleme&L Coinddentally, the spec- 
tral compooents of the ir^f^i^^™ can be formed from the difference of the openings of the 
successive erosioos of tfae image, eq. C7), for algorithm cyrimirarinn, The morphological 
n»onstrnction of the high resolution 3^ data featnres proceeds by selectively addm^ 
^yytf^ skeletal ^< w« f «iwir f eattffes in an interative process from large to small scale based 
on connectivity with the modified Hansdorff metric eq. (12). Fundamentally die skeletal 
feamresareaadedtotbe r econst r ucti o Dbased u pon p r ii icipalsoffnz^ 
ositiao 2.4J2 . Additional r e quir em em s of the morphological ikrWtnn relative to 14-D 
mrvilinfiir feflturw arf addttttH in t**^ foll'^^^g w%^t>wvt* winn nf tfai.< papgf The ftmda- 
mental principal (tf segmenting image structures based on size and connectivity is demon- 
strated in Hgure 3 aiKl Figure 12. 
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ISA Definition : Motphological reconstnictioo based upoo 
fuzzy connectivity , 

R5(0 - (r.iyw) + [s:(g. J) A {r. Rj^.O))] <i5) 
R«4V)«/4. = s;(g./) -0. , 

d«ne<d + E, nEN. for digital images e is 1. 
RJCO is the completed recoostractioa. 

2.6 Morphologkal Scale Spectra of Skeleton Components 

A similarity of the morphological set analysis to Fourier analysis is the ability to represent 
image data with spectral distributioos. The frequency oomem of the znccplK^ 
are however relative to tractable shapes rather than the more abstruse rqvesentation of fre- 
quency forsinusoidalfimctiCHispectraintheFourierdistributioa. The ide^ sh^)e spec- 
tra of images was developed by Matheron ^ to define probabilistic size distribtitioQS of 
shapes for the study of 2-D binary sets on continuous domains. Sena funfaemiore zppMtd 
this tedmique to study the porosity of petrographicandbtok>gical images Tliedevel- 
opment of amotptelogical frequency q>ectrum is based upon the measures of the compo- 
nents of the moiphok)gical skeletal. The frequency magnitude (PSa) of each spect^ point 
(a) of scale or orientation is therefore related to tiie volimie measure of the morphological 
skeletcD oompooents. 

2.6.1 Definition : 

White spectral distribution, 

ps.(«. /) - ^ r./i, - sTo tIIt.^. - IT— /I.] • 

Black spectral distributioQ, 

PS..(fr /> - ^ r./I, - - r./l.) . (17) 

hi, denotes the Lebesque measure on X, \f\i^ J fdm. 

For digital raster formats of pixels or voxels, a is to the integer domain Z and e is 

equal to 1. The spectral disiributioQ of shape for digital images is furthermore ayproximatrd 
with the voltmie o£ the morphological fV^V'^^ tepfesented as the sum of grayscale integer 
voxels intensities : 

2jS^ Ptnrr*****" : Spectral distribotioos for digital data, 

PS.(g. /) - If.>1, - |f,,, fl , white . (18) 

PS..(g. /)- r... /I, - Wi /I, • black , 09) 



ForHgiires3and4tbepatternq>ectraaiialysisof this continix>us biliary inodel of a MRI 
transverse section, ttsough the brain, is ^^''^wpiifi^a qualitatively without high resolution 
features in Hgure 6 for infinitesimal (filatkms and erosions. I^esented in tiiis spectral dis- 
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tribution i$ a positive spectral ccmpcnent where the endrciing extra^oanial ring disappears 
due to erosion at a spectral scale and a negative spectral component wbere the dark ring re p- 
resentationofibeskull disappears due to dilation at anegative spectral scale. For bicsnedical 
image arqp^T^tF"»< or synthetic qxnure radar (S AR) dntter data with rich morphological 
ccmfient the behavior these scale spectral distributions relate a description of the more com- 
plex and &actal nature of biological structures, particularly for the neurological and the vas- 
cularsystems.lliepreservatiooofthisiiatureisanextremelyin^K»tantattri visu- 
al representaticHi and human evaluation of biomedical imagery. 

A morr grnrral Ttitffrw>^arinn of the spectral distribution is related to the Minkowski content 
^ ^ w bTr** can bff gr«ggali»l to ftmrrinng over an nnaye domain with the structuring ele- 
ment base 3, . Ihe result provides a means of measuring £cactal dimension, s, however in a 
more restricted than the Hansdocff dunensioo. (Le. find s soch that log H is bounded) 

2.63 Definition t Minkowski Content |i is a scalar function mapping 
dependent on s € [O^n] of , 
wiiite lower bocmd, 

."7, a^iiiL _ ,.(.) . f./ . / s i . (20) 

black upper bound ^ 

HiAjdk . ^.u) . .To r./- > e <u . (2.) 

2.7 Sixe FQtcrii^ with Morphological Scale Spectra and Skeleton Components 

A further similarity of the morphological size spectral distributions to the Fourier transform 
is the ability to perf onn frequency filtering of the spectral coaqK>nents. Hence it is possible to 
perfonn h=g Mfpff t fT, band-reject, highpass or lowpass filtering of the shape size composition 
of an iixiage. In this inanner a spectnd filter F(a) may be defined over the size spectral ^ 
main« Specifically size fUtering introdoces anew postprocessing methodology for segment- 
ing static tissues from MRA acquisitioos. 

For the S^MRI gnages tbe discriminatin gfeatnreofstatic tissue s of the brain i^ 

propeny across abroad range of neighborhood scales relative to the limited radius of intra* 

f wrnitti art<i4 •fiH vgin.< Filtering tfai^ convcxitv feature over a broad range of scales per- 

mitt the selective eUminadon of static tissues by limiting the scale of morphological recon- 

stnictiOD to within ihb ladins of arteries and veins. In effect this tecfaniq ue is a sel^ 

ture fcal f filtering enhancement technique. To perfonn this preprocessing operation a 

size spectral filter is defined. 

2.7 J Definition : F(a) » 1 for 0 ^ a S v. and F(a) « 0 v a . 

Lt. Morphological rtconstruaion based upon a limited range ofdiametriascale v, vis the 
appnadmme maximum diameter of arteries and vew. 

thf ^ vr^'h fil tf^ thft fTF*^'r^*y p r fwr ti gf of "v^ggf pygra pattem spectrum are selective- 
ly *w*»' J i>^'*^n^f or suppiessioQ to assist visoalizaiioo and image analysis. 
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2.7^ Definition: 

« r./;,.^ + F,(a)s.(g. /) 

wbcc S;(g. /) = F(a)S.(g. /) « 2]F(a)S.(fr /) . (22) 
/ = , is the lim^tPii scale reconstruction. 
%A CorriliDear Feature Segmentation 

The csrviiiDear features of images are discrimi&atBd and enhanced fsxsn. other concave fea- 
tures, such as surfaces and noise, due to their behavicH- relative to 1-D stractnring elements. 
This analysis exploits maximum cooqxarison, by openings and tntntrwai comparison, by clos- 
ing, of 1-D line segments oriented at various angles andscales of kmgth ^'^'^^.Tlte morpho- 
logical skeleton representatico of the image is used due to the reduced dimensianal represen- 
tation of curvilinear forms of the arteries and veins as 1-fD curves of unifbnn diameter. In 
particular a hierarducai rqsresentation by cascading maximmn orientaticHis extends the 
scale of interest for selection of arteiio-venous curvatures. This tedmique may likewise be 
extended to analyze 2-D stxrface fca tuies with 2-D planar sh^ied stmcmring e lemen ts ai 
various angles of orientation. This capability of dimensional fom filtering extends to a 
iw ytfirW py yitr^fvw fnr thf>. npgrarirtfi fif erpftioa^ and complimcntarilv the operation ofdila- 
tion, a c^Mbility of analyzing diminsinnal fonns and rqgesenting data features in the skde- 
ton of ahence lesser dimensional f onn to that of the stracturing elemenL Such as for a spheri- 
cal structuring element the skeleton contains pcnnts, curves, surfaces. Figure 5 lists a chart of 
the discriminated fonns based of the structuring element dimensionality. 

Curvilinear analysis differs frtsn the size and convexity analysis of the preceding segmenta- 
tion schemes primarily due to tlie reduced dimensionality of the stracturing element from 
3-D sphere to a l-D segment and hence the greater sensitivity to the orientations of curvili- 
near strxictnres(Hgixre 8). Ercttion with this stmcturingeleinenteiiinxn noise as 
weU as cnrw and stirfaces points wiiidi do not have parallel tangents to the segment Utili^ 
ing the supremum over all orientations of opening with the 1-D stracturing ele m e n t filters 
noise and the apptpp riate tangent direction of curvilinear features eq. (23). The infi- 
Tnirni ovc aD orientatioQs of closings with the 1-D structuring elenoent dftfinns tlie tangent 
direction. Figure 7, and enhance the curvilinear features of the MRAacquisitio (23). 
The length of the 1-D structuring element is also scaled for improved sensitivity to the cur- 
vatures of curvilinear forms. Figure 8. 

2^1 Definition : The median operations are : 

eretton, H/- V Tf^'f. dilation, F^/- A tT S . 
^eo coeo 

opening, - V /,,^., , closing, - A /^-> . (23) 

<(cd) - denotes the 1 - D oriented strucuring element at - 

angular co within the ^fy^tn Q. 
A deix>tes an infimim over X of a collection of fun ct io ns . 
V denotes an siqscmum over X of a coUectioo of fun cti o ns . 



Analysb mod ScgmcnUtk)n with Fuzzy Operators for ReprcscntationAlsaalization 



Fqnariffn^ of (23) describe modified median filter operations, 
for digital imagery: 

^j- dfi K tf ffi r the 1-D stnicturing element at orientadoa i, 

i € {0 12), denotes the 13 possible digital orientations of a 3-voxel segment within a 

3x3x3 voxel matrix. 

Curvilinear (greduced dimensional analysis furtberoffers ameans of motphological skde- 
tonizatiaQ over axange of scales, relative to the lengths of the varioosly o 
elements. Theanalysis of this skeletondemonsiratescurvatiirei^ curvi* 
ft I f nrm< nf vfftwft mA flitgrieg- Tntnitively. this anal ysis is r ea s flOflh le since smaller carva- 
tures» as qsposed to larger, should be less setisitive to changes in the scale length o^ 
structiirxng element For digital imagery osrvilixsear formis ham the finite dimensions of 
voxels whidi p**mitg cnrvilinear analysis over a broad range of scales. 

2^ DcllnitkMi : Oriented slrrlnnn , 

S;(^,/)«(f;/)\(f;/)^^ . 

Curviliiiear analysis is only perfonned upon the morphc^ogical skeleton, formed fitom eq. 
(7) which exchasivtly contains < 3-D dsneosional fonns. Hiis restrictiGfi upon curvilinear 
analysis is imposed since the varying diameters of arterk^ and vdcns as weU as s^ 
ftf fh#> Krain, in f iw> nri gtnfti MRA data set, are 3-D forms wiiicfa are invariant to the oriented 
1-D structoring elements used for curvilinear analysis (Hgnre S). 

^tfwi fft MtmtMt^ ctari/^rijcgngg Imm the image, Tliermnrrtivitv Criteria isalsouscd to elimi- 
tiafi» PTtT OTMifaii ri jt«ig fgatm»t nf thft motphological skeleton for an analysis focused upon 
1-D curvilinear f eamres. An enhancement of the composite skeleton is thus possible with a 
^i*rh«tq ^T^ eq. (26) similar to the intracranial reconstruction process, eq. (15). 

2^ DcllniiioD: 

CLi;{<,S-) - ( n CLi;.. K,S-)] (26) 

+ s«r.s-) A[rf CLi;..(^,s-)]. 

CLi;{<,S*) - 0. r as r . 
Further eohancement can be provides with median closixig, 

CLiJ(^,s-)- V ([cLi;(e,s-)l)" a?) 

OSrSr 

Fcr 3^ digital MILA Studies equation (26) requires rs: 2 to approximate ar 
vaturesandc« 1 voxel spacing. A simpler approach is to perform only one opening and one 
c k>$ ipfl '?f^* ? ff WH» Tw* for filtfiing ai*** mYLrxrtr-i^m t«pfterivrfy Tliis approach 

however does DOC preserve image featnresweU^i^nr becomes larger doe to the lack ofa 

connecdvity criteria. 



(24) 
(25) 
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S K S*(g, /), is th composite skdexon of the image 

/ dfirived with a spherical stmcturiag elemcsu g. 

of dianietxic scale S the ip?^™mn ctiamftfT of arteries and veins. 

2^ Definition: 

CL2(^S-) « S* A (S;,f , where S* « S'Cg, /) . (28) 

T?/ |narinn (2?) yields noisc suppression and enhancement of the limitrd scale moxphoiogical 
cwwri typicaUy r » 1 or 2 fcff 3^ digital MRA studies. In general to retrieve the true 
A\mtm eL tyjpfMgntatiop of tfac arteries and veins ti>e results of eauadons (26. 27) are used to 
reconstruct a filtered represeatation of the image, 

2AS Definition: 

r;(/) - [ n r;*.(/)J + s; A cLn(^s-) . s; « s; (g. /) . (29) 

R2(/) is the iwni»i*rf scale reconstroction f or r » 0 . 

3.0 Remits 

n*n ycfta v> mrtr phftif y^i anal yse nf a !qmi pmcno image scction tfarongh the brain was used 
to illusaatethe elementary aspects of morphological description. InHgure lOan example of 
segmenting the intracranial brain tissue from the extracranial tissues of tiie scalp is pres- 
ented. Figure 10c in comparisoo with Figure 10b. The method utilized is similar to die ideal- 
ized version shown in Figure 3, however* the fuzzy connectivity reconstruction tnThni que. 
eq. (15), was taed to restore sinall scale inoq)hological features of the brain after size filter- 
ing. This same technique was extended to 3 dimensions to segment intracranial features for 
tmobscured sor&oe visualization of the brain. Figure 12. 

T5> ^V>wK «nT*ril*^ <|«y^hrtlng^/>al «W1<»tri*^ anrf itg hghavifrr thft xkelfctal mmponents S^^ of the 

image were derived witheq. (7) tip to themaxhntmi diameterd which eroded the image to a 
m m o^CT>e v^ve TNm pisel simcmring elements were tised in a ra vartwi su cressi on of 
erosions togeneratethespectnmiiof iteoessary scales, eq. (4). The shapes of the stmctiiring 
elements were sdected to ap pr oxim ate a drde, one of a diamond shape and the other of an 

nr^fl ^yl ^^i p#^ ^ cHrttilrf f>fwpHa«-M>^ thnr thf> fMm.<migrinn nf the image from its mof* 

phological i^^^^rw^ components, eq. (10), is lossless relative to the original image. The sum 
of th ffi ft skeletal coaqxneats, eq . (8), are visible in Rgure 1 0^ This Hgure illustrates a 1-D 
curvilinear fdm of the composite which is a reduced representation of the 2^D 

moqAology of the imaged brain section. Furthermore the integral volume of the skeletal 
coaqwoentseq. (18), over a sp e ui urn of diamrtm prodDces a size spectral distribtition. Fig* 
urell, Mlirh itf TP^**" ^^ ^ Hgm<rfnr fr nr franinnal dimgnsinn of tfac brain section relative 
to a logarittonic •^^^ In particular the first slope over a range of the smallest structuring 
oUwvnt diametm tp<fi r? i y'f a 1 .65 Hausdocff dimmtinn wlxich is basic to the morphology 
the brain gyxirepKCsentBd in this 2p-D section. By selecting a particular setof**boi^**,Bg- 
urelOd,fiom the skeleton it is possible to independently reoonstnia the image features of 
the bones. Figure lC)e, or to see the image with the feamres of the bones iiussing,Figt2re lOf. 
Thisexample partini^^^y tgpgjtgtit^ tiig muMral and fazzy Teprrsuntaiinn d image features 
with grayscale morphology. The techniqoe introdnced for the analysis of the morphological 
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^vrii>^r« of MRA acquisitions, cq. (26, 27, 28), provides ameans of selecting these skeleton 
bones based f eatores of one dimensiooal curvilinear forms. These broad c^>alniities of the 
moq^bological fV-^^^^nn to represent imagery without distortion, inserve connectivity, iso- 
late diametric featmes, and to fiirthcnnofc represcm the natural or 
posed the mcvphological skeleton as an ideal and extremely attractive representation for ar- 
teries and veins, finun MRA acquisitions. 

Fignxe 13 demoiistratesainagnetic resonance angiogis^y study of the caroddarteriesof the 
yti»Hr Thf yn'^ngy an* selfctp^ rfflaf^a tn flow velocity directions toward the extremity of the 
ft fifF* hmiy An mdesirable artifact from this imaging tcdmiqne is that arterial flow 
obliqpf to tfag irlfCtnT dtregtinn isditrnrtfri resulting in a staircase e£FecL The objective 
t»ith gfnHy u tn jgihance the arteries within the 3-D volnme of the data. The i madnhrr - 
aN>H rfata w mtttan y nh«grved with a transparency projection algorithm, performed through a 
succesriooofpnvection angles. TtefirstinetlKxi to exploit the curvilinear features of arter- 
ies and veins is to lepresent arterio-venoos features in a unifonn diametric rqnesentation. 
This representation is rendered by the moqibological skrlrffmirarinn o£ arteries and veins 
whidi definesacentialaxis of the arteries Le..Fxgurel3b.Tl>e8rteries of the skeleton repre- 
sentation arc approximately of a single voxel diametric width to snnplify cunolinear analy- 
sis of the arterio-venous features. Hiis single voxelrepfesentation allows features of curva- 
ture and length lobe axialyzed independent of die dianietric scale of the artexy or vein. The 
limi ted moq^ogical reconstmctioD, eq. (22), was used to enhance this smdy by elim- 
inating static tissue over a broad range of scales followed by a^ilication of eq. (26, 28) to 
f^ttntnfl ty c^oUary feanses. Figure 14c versus Hgure 14f. Of particular importance it the 
fact ttiat the diameten of the arteries and vetiis are unchanged by this scale selective tech- 
nique . AdditkmaUy, the **stair case** artifact has been eliminated firom di^ 
al aoqnisitioos. 

The visnalizatians of Figure 14 were acqmredatacocrserresolution than Figure 13, howev- 
^ t*"^ T n! "* mrv^^^ thfthmadgTTange of the carotid arterial tree wiiich feeds thebaseof 
th^i^tn Ttw> Kami artmaifeaniM of ttie brain were nihanr^ inFigtce 14a, by the elimi- 
nation ofposterior and anterior coronal slices for a region of interest (ROD smdy. To elimi- 
nate the e xp er t intensive and cooqiuter interactive trrhnignrs of ROI selectioiv die auto- 
Tnyt'^H^t^^ tfamfaoftbiypflr^'^*^*^^^ 

it.w^«-»>^ »i«i^irrarir>n nf tift^rifli tiM Ongg ayaineo. (22^ is rliminatrs image features of 
Hiqmii»»w igwyrthflti thg ftrteriga flnd veins, thus onlv alimitcd scale moTPholog- 

i/^i gv^wrw* t er mfifflnrti™ gurvilinear segmentation of the limited scale mor- 

ph^lrfcgjya l •V^l*'*^ «^!** ^ ^ •^^^ j>ltfn W>ati>c OJ^ pnint iiftt.<g and limits artgrial car. 

vamtes with fe*^ segments of various lengths and orientations. Thus segments of 3 and 5 
voxels in length, or radios, r « 1 and 2, were used with values of oriextfation in 1 3 

directk»s» eq« (23). This segmented tVr^^ then provides a basis to morphologically re- 
mwjtmwrt ^ ffthatiggd data vnlmng fortTeprojection and vitnaliTation. Hgures 14b aiKl 14e 
show f» ^ w *p ^ff of the filtered mmpbological ikrtrtm axinometrically projected along 2 dif- 
ferent angles of direction. Figures 14c and 14f are the projections of the filtered MRA vol- 
ume. Note the preservation of snudl diametric features between prorjeoions of the original 
volimie 11a and the filtered vohnne 11c, and note the radical ehminatioo^f the obscuring 
capillaries from Hgure 14d demonstrated in Figure 14f. 

Figure IS demoostxates a smdy for visualizing the arteries and veins of the bead. Againa 
dtrectiooal fkvw field was selected to enhance blood flows through tike transverse pLa^ 
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wardsthesupenorof the bead. Since aneries and veins have no particular dro in the 
head both types of flow are realized by this imaging method. The elimination of extracranial 
veins as well as c^nllaries is ttie objective for this smdy, particularly since intracranial arter- 
ies axKi veins are of paramount importance lot plaxming neurological surgeries. To farilitarr 
this objective, eq. (1 5) has shown the capability of segmenting intracranial features for unob- 
scured 3-D visualization of the brain. Figure 12, and now may be used to include the seg- 
mentatioa arteries and veiiis within the proxirriity of the static brain tissues. Rgures ISa and 
ISb demonstrate two projected directions of the MRA smdy, extracranial c:q}illaries obscure 
f vt^^^ltTati rm nf ftw mtraeranial arteries and veins and must he e lhnmatf/l Equarion (1 S\ 
was used to seginent the static tissues of the brain based iipcQ a size and ocxmectivity ch 
which pennitted the tmobscured visualizatioD of brain anatomy, Hgme 14, Figure ISeX 
The static tissue representation defined with this xnethod was tiien tised to xiiask t^ 
Dial arteries and vdns forftmhercnrvilinearmhanrfmept eg. (15) and limitrri scale filtering 
eq. (22), and eq. (26, 27) whicfa resulted in Hgure 12cd. Comparing Figure 15a,b with Hg- 
ure IScd. we observe that the intracranial arteries and veins are segmented successfully 
without the overiying extracranial features. 

In Figure 16 a ^w^pw scene is analyzed consisting of foliated trees, tree stump, road and 
telq]^ooe poles. The origiiial iiziage is demcnstraiedin the iipi^ 

structnring ekment is used, operator g in eq. 21, to devdqp a morphological skeletoiu and 
then iccoiistrtia the iiiiage based cm fuzzy conricrtivity of the nkri^ 
note the preservation of foliage details and eliiiunatiOD of inanmade features such as the ro^ 
cign anrf lani^ *trip m th^ mftri In thft Inwer left a 3x3 median opening operation - consistinp 
of 4 1^ oriented structuring ^^b-^twmt^, was exrq)loyed in a succession of 20 rasradrd opera- 
tions This operation preserved the l-O features of the image, including man-made fea- 
tures of ttie sign and strip on the road, aiul eliminated the f eanses of leaves from the 
trees. The leaves have a cfaaraaeri Stic dimcMinn less than one in the im 
by the 1-D oriented operations of the median erosion, eq. (23). In thelowerri^t a median 
dilation, eq. (23), consisting of four 1-D oriented structuring elements, was used on the 
eroded image to fill the textual features segmented by median erosion. This mrdian opera- 
tion of erosicDS followed by dilations rcstilts in a more convex representation of image fea- 
tures such as leaves. 

In Hgure 17, a synthetic aperture radar image is analyzed coosistingof several specular tar- 
gets in a spedde ixaised Gutter background. The CBighial linage is deiricxist^^ 
left Successively morphological rf%^\m dilation cyeraticais are used to fill the fractal fea- 
tures %A the dut&er and leave the impulsive ^^^^^^ details of the man-made targets unal- 
tered « 

4*0 CoociBsloD 

Matheinaticalmoq>hology has provided a powerful segzzien^on tool for the unobscured 
presentation of 3-D brain anatomy ^ ^ and now can ^litnmjtt^ noise and static tissues for a 
clear demoostratiop of the intracranial arterial tree fttm MRA ac q nw 
fiw structures of arteries and veins are preserved by the curvilinear feature segmentation 
with mathematical morphology. 

A theoretical development of a battery of morphological analysis t r rh niqiirs have beende- 
veloped and applied to biomedical and remote sensing image data. These trchni qnm anal yti- 
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cally capture the dimensionality, fractal natare and fusy set structme of features from the 
image acquisidons. A variety of structming elements, digitally familiar to spheres, planes, 
and lines i^^"^ over a 3-D dcmain, have been utilized over a spectrum of scales and 
o riitinirim g tn cggment featnres of gmilar or complementary dimensions and shape. In par- 

^1 fiilar, th e rfmip ]^>mgntary Hitw^n wimil fg^tmift.^ nf fractftU am p <ie M ^ved ft.^ deviatjons in the 

moipboiogical skeleton across a wide range of scales relative to the riiamrrrr of the convex 
stractoring elements. 

The image analysis of mathematical moqdiology has demonstrated an ability to segment 
3 D biomgdical tmngr framm f?f ma£Tv!^^ r^*^™^ tmagi*fy orgmgitinng atyf mnm\ 
c hy fti>r f eanses from remocely sensed S AR images for an tmobscored visualization of fea- 
tures of interest £romaxinainetTicsinn projections. In comparison with i n tensity histogram 
and thresholding techniques the distribudoos of sh^ spectra, eq. (16, 17), of data are 
!hmh <>»^ ma pp ed as anew analysis and visualization tool, e.g. eq. (22), just as intensi- 
ty histograms are thresholded and msqsped. 

Three pr o gres sively improving methods for segmentation of die curvilinear features of ar- 
teries and veixis from MRA acquisitions were XDtroduced. A successive ^li^ 
methodsiliustratBdanevolutionaryrefixKment of algorithms to segmentintracranial arteries 
and veins ton MRA smdies. One of these methods alternately jsovided a means of seg- 
menting brain tissues £ram the involvement ctf extracranial tissm eq. (IS). Alternately, a 
technique established a means of siq^xessing static tissue features £ram MRA projections 
eq. (22), and enhancing curvilmear features eq. (26, 27, 28). This technique successfiilly ad- 
dressed the segmentation requirements of the peculiar curvilinear nature of MRA acquisi- 
tions, a resplendent example is the comparison ofHgures 14dandl4f,andatestamentofits 
accuracy is exenqilified in the preservation and enhancement of delicate curvilinear fea- 
tures, comparison of Hgure 14a and 14c. 

The premise of imagery as a representation over a Borel field topology in mathematical 
morphology f pr*i«t<i>»* the analysis of imagery with fuzzy set al g ebra and offers a catalyst for 
fiiitheradvanoernentintothefie]dofdatahiterp(^onarKlreprojection.I^ MRI 

arqtitct»tr»%< ttfiH S AP **"ftfl^ ^ a wpr w^wirarirm ftf natitral nrgan«m» with a p r npert y of 

selfdmilatityat various scales ^whidi is consistent to a fractal description ^ The utilization 
of the algebra of mathmiariral inorphology makes possible die direct analysis and seg- 
mentation of fractal image features relative to spectral distributions of size and shape from 
thenxspbological dceleton^ as presentedinttiis paper, eq. (16, 17),Hgiares 10,11. Thisca9>a- 

bility of fractal ir p f »taKH<h^ o Ka«« fnr th<» moct Hteral and mfofmari ve preaeiw 

tation and segmentation of biomedical imagery. 
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figure? 
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Hgm 9 : Ventricles and blood sinuses of the bndn. Low contrast ventricular 
system is demonstraied tarn overlying subdmal cerebral-spinal fluid (CSF) 
by enbancementof dilation. Blood sinnses, inner brain taematnma , andhemor- 
xliaging superior to tbe brain axe demonstrated by enhancement of erosicm. 
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(a) (b) (c) 




(d) (e) (f) 



Figoxe 10 : An anay of image processing operations on aprocon spin acquisition 
of a brain secticm (a) Morphological skeleton with a selected enhancement of 
a 'lx>ne** (b) Segmented neurological tissue (c) Original MRI scan section (d) 
Selected feature bone of the skeleton (e) Morphologically reconstructed gyri 
from selected **bCHie** of skeleton (f) Morphologically reconstructed brain tissue 
without gyri bone feature from skeleton* Note nattval continuity of texture on 
the boundaries wliere features are removed. 
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unit slope 



log(voxel diameter) 



Figme 11 : Log-log distribmicn of size ^lecaa relative to the volume measure 
of each ^>ecgal ccmqxmgnt of the morphological sk-rlr ton. In the smallest scale 
region of the distribution the slope corresponds to the fractal dimension of the 
gyri of the brain scan section of Figure 10. Utilizing the liist eight voxel diame* 
tenforaleastmeanssquare qnadradic fit estimates the gyri to have an 1.65 di- 
mensional repceseatatira within the 2-D scan section. 
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Hgnte 12 : Sum projectioos (SP) of a S-D morpbologically segmented and xe* 
constructed MRI brain scan. Brain surface is grayscale gradient f Tihanml and 
shaded with a single light soiffce at a direction fma infinity, (a) Lateral sagittal 
SP (b) Nearpostenca-cotonal SP (c) InfeikH' transverse SP (d) Superior trans- 
verse SP. The method of visualization is based oa a voxel based stirface en- 
hancement ^'^Ml>e classification of internal brain anatomy has also been iq>- 
pxoached with multi-modal fusion techniques ^ , 
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(a) (b) 




(c) 



Figure 13 : Coronal sum projections (SP) of variously 3-D processed MRA ac- 
quisition of the carotid arteries, phase samrated, (a) Surface enhanced SP(b) 
SP of the morphological skeleton demonstrating curvilinear features, note the 
uniform diametric representation (c) SP of the 3-D acquisition. 
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(a) 



(b) 



(c) 



1^ 



(d) 



(e) 



(f) 



Figure 14 : Demonstration of variously 3-D processed MRA acquisitions with 
ccMTonal sum projection (SP). top row is coxonal sum projection (SP), bottom 
row is near lateral SP (a, d)SP of original data, posterior and anterior sections 
removed to reveal basal arteries^ note the obscuring facial ca pill a ri es in (d)> (b, 

e) SP of curvilinear (CV) *^hyir^ limited morphological slceleton (LMS). (c 

f) SP of mcnpbologically reconstructed data from CVE limited morphological 
skeleton, note the preservation of rf'*^«rt>*T arterial features and the filtering of 
facial capillaries. 
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(a) 



(b) 



(c) 



(d) 



(e) 



(f) 



Hgue IS : Sam projecdcms (SP) of blood vascular and static brain tissue fea- 
tures, left column is an anterior coronal SR light row is a near left lateral SP (a, 
b) SP of MRA acquisition with limited scale reconstruction, (c, d) SP of CV 
enhancement of intracranial blood vasculamre. Segmentation of intracranial 
vasculature was realized with segmented static brain tissues, (e, f) Fused surface 
shaded visnalizadon <tf vasculature and stadc brain tissues. 
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(a) (b) 




(c) (d) 



Figure 16 : Median filtenng on 2f-D imagery provides a means of removing 
fractal rfTTfi*^ratjii f eamres and **filling** textures for an alternate texture map- 
ping. Tbc original image (a) is deccmposed into mrdian erosico skeleton com- 
ponents and reconstructed with eq. (26) (b). In (c) the original image is pro- 
cessed with 20 ra sr^''^ operations of mr^^™ erosion. In (d) the image of (c) 
is processed with 20 median dosing operations. 
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(a) 



(b) 



(c) 



(d) 



Hguxe 17 : Median faltering on 2*D imagery provides a means of removing 
fractal dimesianal featnres and **filling** textures for an alternate texture m^ 
ping. Tbe original synthetic 2pmm radar (S AR) image (a) is median closed 
ivith 8 opetations in (b). 32 operations in (c) 128 operations in (d). 
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WHAT IS CLAIMED IS: 

1. A method for isolating anatomical structures 
contained within a three-dimensional data set, the method 
comprising the steps of: 

5 a) forming a morphological skeleton of the 

three-dimensional data set; 

b) selecting a seed data point within the 
morphological skeleton, the seed data point being 
contained within the desired anatomical structure; and 
10 c) utilizing fuzzy connectivity to define 

additional data points of the desired anatomical 
structure so as to reconstruct substantially only the 
desired anatomical structure; 

d) wherein reconstruction of substantially only 
15 the desired anatomical structure facilitates viewing 

and analysis thereof. 

2. The method as recited in Claim 1 wherein the step 
of forming a morphological skeleton comprises recursive 
opening and erosion of the three-dimensional data set so as 

20 to form a plurality of residuals which define the 
morphological skeleton. 

3. The method as recited in Claim 1 wherein the step 
of forming a morphological skeleton comprises utilizing a 
generally spherical structuring element in recursive 

25 opening and erosion of the three-dimensional data set. 

4. The method as recited in Claim 1 wherein the step 
of selecting a seed data point comprises positioning a 
cursor at a desired point on an image being displayed on a 
monitor . 

30 5. The method as recited in Claim 1 wherein the step 

of using fuzzy connectivity to define additional data 
points of the desired anatomical structure comprises 
defining connectivity based upon the size and shape of a 
structuring element utilizing a fuzzy generalization of 

35 mathematically defined distances between sets of data 
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points as a criterion* 

6. The method as recited in Claim 1 wherein the step 
of reconstructing substantially only the desired anatomical 
structure comprises recursive dilation and closing of a 

5 selected portion of the morphological skeleton. 

7. The method as recited in Claim 1 wherein the step 
of utilizing fuzzy connectivity to define additional data 
points comprises defining connectivity based upon the use 
of a generally spherical structuring element for defining 

10 distances between adjacent data points. 

8- The method as recited in Claim 1 wherein the 
three-dimensional data set comprises a data set generated 
by a device selected from the group consisting of: 

a) a magnetic resonance imaging device; 
15 b) a computer aid tomography device; and 

c) a positron emission tomography. 

9. A method for facilitating analysis of three- 
dimensional medical images, the method comprising the steps 
of: 

20 a) separating anatomical features from one 

another via dimensional filtering; and 

b) storing desired separated features for 
analysis . 

10. A method for facilitating analysis of three- 
25 dimensional medical images, the method comprising the steps 

of: 

a) defining anatomical features having less 
than one dimension from the image via first 
dimensional filter and storing those features; 
30 b) defining anatomical features having less 

than two dimensions from the image via a second 
dimensional filter and storing those features; 

c) defining anatomical features having less 
than three dimensions from the image via a third 

35 dimensional filter and storing those features; and 
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d) storing remaining features having three 
dimensions in a fourth memory; 

e) wherein anatomical features of different 
dimensionalities are thus separated from one another 

5 so as to facilitate isolation of desired anatomical 

structures. 

11. A method for reconstructing an anatomical 
structure from a morphological skeleton, the method 
comprising: 

10 a) selecting a seed data point within the 

morphological skeleton, the seed data point being 
contained within the desired anatomical structure; and 

b) utilizing fuzzy connectivity to define 
additional data points of the desired anatomical 

15 structure so as to reconstruct substantially only the 

desired anatomical structure; 

c) wherein the use of fuzzy connectivity 
results in reconstruction of substantially only the 
desired anatomical structure and substantially lacks 

20 surrounding tissue. 
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